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' Abstract 

Ci 

Weak gravitational lensing is responsible for the shearing and magnification of the im- 
^ , ages of high-redshift sources due to the presence of intervening matter. The distortions are 

^ I due to fluctuations in the gravitational potential, and are directly related to the distribution 

of matter and to the geometry and dynamics of the Universe. As a consequence, weak grav- 
itational lensing offers unique possibilities for probing the Dark Matter and Dark Energy 
in the Universe. In this review, we summarise the theoretical and observational state of the 
subject, focussing on the statistical aspects of weak lensing, and consider the prospects for 
weak lensing surveys in the future. 

Weak gravitational lensing surveys are complementary to both galaxy surveys and cos- 
mic microwave background (CMB) observations as they probe the unbiased non-linear mat- 
ter power spectrum at modest redshifts. Most of the cosmological parameters are accurately 
estimated from CMB and large-scale galaxy surveys, so the focus of attention is shifting to 
understanding the nature of Dark Matter and Dark Energy. On the theoretical side, recent 
advances in the use of 3D information of the sources from photometric redshifts promise 
greater statistical power, and these are further enhanced by the use of statistics beyond 
two-point quantities such as the power spectrum. The use of 3D information also alleviates 
difficulties arising from physical effects such as the intrinsic alignment of galaxies, which 
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can mimic weak lensing to some extent. On the observational side, in the next few years 
weak lensing surveys such as CFHTLS, VST- KIDS and Pan-STARRS, and the planned 
Dark Energy Survey, will provide the first weak lensing surveys covering very large sky 
areas and depth. In the long run even more ambitious programmes such as DUNE, the Su- 
pernova Anisotropy Probe (SNAP) and Large-aperture Synoptic Survey Telescope (LSST) 
are planned. Weak lensing of diffuse components such as the CMB and 21cm emission can 
also provide valuable cosmological information. Finally, we consider the prospects for joint 
analysis with other probes, such as (I) the CMB to probe background cosmology (2) galaxy 
surveys to probe large-scale bias and (3) Sunyaev-Zeldovich surveys to study small-scale 
baryonic physics, and consider the lensing effect on cosmological supernova observations. 

Key words: Gravitational Lensing 
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1 Introduction and notations 



Gravitational lensing refers to the deflection of light rays from distant sources by the grav- 
itational force arising from massive bodies present along the line of sight. Such an effect 
was already raised by Newton in 1704 and computed by Cavendish around 1784. As is 
well known, General Relativity put lensing on a firm theoretical footing, and yields twice 
the Newtonian value for the deflection angle (ItiI) . The agreement of this prediction with 
the deflection of light from distant stars by the Sun measured during the solar eclipse of 
1919 (70) was a great success for Einstein's theory and brought General Relativity to the 
general attention. The eclipse was necessary to allow one to detect stars with a line of sight 
which comes close to the Sun. 



In a similar fashion, light rays emitted by a distant galaxy are deflected by the matter dis- 
tribution along the line of sight toward the observer. This creates a distortion of the image 
of this galaxy, which is both sheared and amplified (or attenuated). It is possible to dis- 
tinguish two fields of study which make use of these gravitational lensing effects. First, 
strong-lensing studies correspond to strongly non-linear perturbations (which can lead to 
multiple images of distant objects) produced by highly non-linear massive objects (e.g. 
clusters of galaxies). In this case, the analysis of the distortion of the images of background 
sources can be used to extract some information on the properties of the well-identified 
foreground lens (e.g. its mass). Second, cosmic shear, or weak gravitational lensing not as- 
sociated with a particular intervening lens, corresponds to the small distortion (of the order 
of 1%) of the images of distant galaxies by all density fluctuations along typical lines of 
sight. Then, one does not use gravitational lensing to obtain the characteristics of a single 
massive object but tries to derive the statistical properties of the density field as well as the 
geometrical properties of the Universe (as described by the cosmological parameters, such 
as the mean density or the curvature). To this order, one computes the mean shear over a 
rather large region on the sky (a few arcmin^ or more) from the ellipticities of many galax- 
ies (one hundred or more). Indeed, since galaxies are not exactly spherical one needs to 
average over many galaxies and cross -correlate their observed ellipticity in order to extract 
a meaningful signal. Putting together many such observations one obtains a large survey 
(a few to many thousands of square degrees) which may have an intricate geometry (as 
observational constraints may produce many holes). Then, by performing various statisti- 
cal measures one can derive from such observations some constraints on the cosmological 
parameters as well as on the statistical properties of the density field over scales be tween a 



^sity nei( 

few arcmin to one degree, see for instance Refs. (1 1 971 : 1 1 901 : [l7t I2191 : 12991 : 1 1 9 11 : 12441) . 



Traditionally, the study of large scale structures has been done by analyzing galaxy cata- 
logues. However, this method is plagued by the problem of the galaxy bias (i.e. the dis- 
tribution of light may not exactly follow the distribution of mass). The advantage of weak 
lensing is its ability to probe directly the matter distribution, through the gravitational po- 
tential, which is much more easily related to theory. In this way, one does not need to 
involve less well-understood processes like galaxy or star formation. 
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In the last few years many studies h ave managed to detect cosmological shea r in random 
patches of the sky 0; S; S; M, [112; EH; H; H H; llisl;!^; S H liooh . 



While early studies were primarily concerned with the detection of a non-zero weak lensing 
signal, present weak lensing studies are already putting constraints on cosmological param- 
eters such as the matter density parameter and the amplitude Og of the power- spectrum 
of matter density fluctuations. These works also help to lift parameter degeneracies when 
used along with other cosmological probes such as Cosmic Microwave Background (CMB) 
observations. In combination with galaxy redshift surveys they can be used to study the 
bias associated with various galaxies which will be useful for galaxy formation scenarios 
thereby providing much needed clues to the galaxy formation processes. For cosmologi- 
cal purposes, perhaps most exciting is the possibility that weak lensing will determine the 
properties of the dominant contributor to the Universe's energy budget: Dark Energy. In- 
deed, the recent acceleration of the Universe detected from the magnitude-redshift relation 
of supernovae (SNela) occurs at too late redshifts to be probed by the CMB fluctuations. 
On the other hand, weak lensing surveys offer a detailed probe of the dynamics of the Uni- 
verse at low redshifts z < 3. Thus weak lensing is among the best independent techniques 
to confirm this acceleration and to analyze in greater details the equation of state of this 
dark energy component which may open a window on new physics beyond the standard 
model (such as extra dimensions). 

In this review we describe the recent progress that has been made and various prospects 
of future weak lensing surveys. We first describe in Section [2] the basic elements of the 
deflection of light rays by gravity and the various observables associated with cosmolog- 
ical weak gravitational lensing. In Section [3] we review the 2-point statistics of these ob- 
servables (power- spectra and 2-point correlations) and the problem of mass reconstruction 
from observed shear maps. Next, we explain in Section |4] how the knowledge of the red- 
shift of background sources can be used to improve constraints on theoretical cosmological 
models or to perform fully 3-dimensional analysis (3D weak lensing). Then, we describe in 
Section|5]how to extract further information from weak lensing surveys by studying higher- 
order correlations which can tighten the constraints on cosmological parameters or provide 
some information on non-Gaussianities associated with non-linear dynamics or primordial 
physics. We turn to the determination of weak lensing shear maps from actual observations 
of galaxy images and to the correction techniques which have been devised to this order in 
Section[6l In Section|7]we discuss the numerical simulations which are essential to compare 
theoretical predictions with observational data. We describe in Section [8] how weak lens- 
ing surveys can also be performed at other wavelengths than the common optical range, 
using for instance the 21cm emission of first generation protogalaxies as distant sources. 
We present in greater detail the weak lensing distortion of the CMB radiation in Section |9l 
In Section [To] we also discuss how weak lensing can be combined or cross-correlated with 
other data sets, such as the CMB or galaxy surveys, to help constrain cosmological models 
or derive some information on the matter distribution (e.g. mass-to-light relationships). Fi- 
nally, we conclude in Section [TT] To help the reader, we also give in Tables [I]-[3]below our 
notations for most coordinate systems and variables used in this review. 
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Table 1 

Notation for cosmological variables 



total matter density in units of critical density 



reduced cosmological constant Ha 

reduced dark energy density Cl^e 

mean comoving density of the Universe p 

Hubble constant at present time Hq 

Hubble constant at present time in units h 
of lOOkm.s^'.Mpc^i 

rms linear density contrast in a sphere Og 
of radius Sh~^ Mpc 



Table 2 

Notation for coordinates 



metric ds^ = c^dt^ - a^{t)[dx^ + ©^(de^ + sin^edcp^)] 

speed of light c 

scale factor a 

comoving radial coordinate %, r 

comoving angular diameter distance © 

comoving position in 3D real space x, r, 2) 9) 

comoving wavenumber in 3D Fourier space k, {k\\,k±),{k\\,£/'D) 

bend angle a 

deflection angle 69 

image position on the sky 9 

flat-sky angle (9i,92) 

2D angular wavenumber £, {£x,'^y) , (^i ) h) 
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Table 3 

Notation for fields and weak-lensing variables 

gravitational potential <J> 

lensing potential (|) 

shear matrix ^ 
amplification matrix 

weak-lensing convergence K 

complex weak-lensing shear Y = Yi + ^ Y2 

shear pseudo- vector Y = Yi + li^y 

tangential component of shear YnY+ 

cross component of shear Yx 

weak-lensing magnification /i 

angular filter radius 6^ 

smoothed convergence, smoothed shear ic, y 

weak-lensing aperture mass Map 

3D matter density power spectrum P{k) 

2D convergence power specttum ^k(^) 

2D shear power spectrum Py{^) 

two-point correlation ^ 

3D density contrast bispectrum B{ki,k2,kz) 

2D convergence bispectrum Sk(^i,^2,^3) 

probabiUty distribution function of the smoothed convergence S'kCk) 



9 



2 Weak Lensing Theory 



2.1 Deflection of light rays 



We briefly describe here the basic idea behind weak gravitational lensing as we present a 
simple heuristic derivation of the first-order result for the deflection of light rays by gr avity . 



For a rigorou s derivation using General Relativity the reader can consult references (I152I : 



17l: l232l : l252h 



We assume in the following that deflections angles are small so that we only consider 
first-order terms. This is sufficient for most applications of weak lensing since by defini- 
tion the latter corresponds to the case of small perturbations of light rays by the large- 
scale structures of the universe. Let us consider within Newtonian theory the deflection 
of a photon with velocity v that passes through a small region of space where the gravita- 
tional potential 4> is non-zero. The acceleration perpendicular to the unperturbed trajectory, 
= — V^4>, yields a small transverse velocity = — / d?V^^>. This gives a deflection 
angle 6c = v^/c = — / d/Vj^^>/c^ for a constant velocity |v| = c. As is well-known, General 
Relativity simply yields this Newtonian result multiplied by a factor two. This deflection 
changes the observed position on the sky of the radiation source by a small angle 50. For 
an extended source (e.g. a galaxy) this also leads to both a magnification and a shear of 
the image of the source from which one can extract some information on the gravitational 
potential 4>. If the deflection takes place within a small distance it can be taken as instan- 
taneous which corresponds to the thin lens approximation (as in geometrical optics) as dis- 
played in Fig. [U Besides, in cosmology transverse distances are related to angles through 
the comoving angular diameter distance T) given by: 



csinjf (|i-an,-aA|i/2//ox/c 



where sin^- means the hyperbolic sine, sinh, if (1 — Q.^^ — ^a) > 0, or sine if (1 — — 
Q.a) < 0; if (1 - Q.jn - ^a) = 0, then © (%) = % (case of a flat Universe). The radial co- 
moving distance % measured by a light ray which travels from a source at redshift z to the 
observer at z = is given by: 



y = ^r , (22) 

HoJo v/aA+(i-a„,-fiA)(i+z02+f^m(i+zO^' 

where z' is the redshift along the line of sight. Note that % measures both a spatial coordinate 
distance and a travel time. Here we also introduced the Hubble constant Hq and the cos- 
mological parameters Q.^ (matter density parameter) and Q.a (dark energy in the form of 
a cosmological constant). Therefore, the source appears to have moved in the source plane 
over a comoving distance © (%i')50 = —'D{Xs — X)oc as can be seen from Fig. [H where d 
and 60 are 2D vectors in the plane perpendicular to the unperturbed light ray. Summing 
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source plane 

lens plane 

observer plane i 1 




Fig. 1 . Deflection of light rays from a distant source at comoving radial distance %s by a gravita- 
tional potential fluctuation at distance %. For a thin lens the deflection by the angle a is taken as 
instantaneous. This changes the observed position of the source by the angle 66, from the intrinsic 
source direction 6^ to the image direction 6/ on the sky. 

up the deflections arising from all potential gradients between the observer and the source 
gives the total shift on the sky: 

50 = 0,-0, = -/ dx )^',^' V±^{X). (2.3) 



where 0, is the intrinsic position of the source on the sky and 0/ is the observed position. 
However, generally we do not know the true position of the source but only the position of 
the observed image. Thus the observable quantities are not the displacements 60 themselves 
but the distortions induced by th ese defl ections. They are given at lowest order by the 



symmetric shear matrix (|146Hl7Hl37() which we define as: 



- 30,, - c2 io 'DiXs) ^^•''^ 

Eg. (12. 41) follows from eq. (|2.3l) if we note that a change of angle d0 for the unperturbed 
light ray corresponds to a transverse distance D (%)d0 in the lens plane where the gravita- 
tional potential ^» produces the gravitational lensing. The reasoning presented above clearly 
shows that Eq. (|2.4l) uses the weak lensing approximation; the derivatives V/Vj^>(%) of the 
gravitational potential are computed along the unperturbed trajectory of the photon. This 
assumes that the components of the shear tensor are small but the density fluctuations 5 
can be large (146). We can also express the shear matrix in terms of a lensing potential 
(|)(0;%i.) (also called the deflection potential) as: 

^'.j = tu with mx.) = I I^X ^(%,^(X)9). (2.5) 

The expression (12. 5|) is formally divergent because of the term 1/© (%) near % = 0, but this 
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only affects the monopole term which does not contribute to the shear matrix (indeed 
derivatives with respect to angles yield powers of © as in Eq. (l2.4l) ). Therefore, we may set 
the constant term to zero so that (|)(0) is well defined. Eg. (12. 51 ) clearly shows how the weak 
lensing distortions are related to the gravitational potential projected onto the sky and can 
be fully described at this order by the 2D lensing potential (|)(0) . Thus, in this approximation 
lensing by the 3D matter distribution from the observer to the redshift Zs of the source plane 
is equivalent to a thin lens plane with the same deflection potential (|)(9). However, from the 
dependence of <^{B;Xs) on the redshift Zs of the source plane we can recover the 3D matter 
distribution as discussed below in Section |4l Note that weak lensing effects grow with the 
redshift of the source as the line of sight is more extended. However, since distant galaxies 
are fainter and more difficult to observe weak lensing surveys mainly probe redshifts ~ 1 . 
On the other hand, this range of redshifts of order unity is of great interest to probe the dark 
energy component of the Universe. Next, one can also introduce the amplification matrix 
of image flux densities which is simply given by the ratio of image areas, that is by the 
Jacobian: 

^ = |i=(5,+^,)-^=(' - - I, (2.6) 




39/ 

which defines the convergence K and the complex shear y = yi + /Y2. At linear order the 
convergence gives the magnification of the source as /j = [det(.;^)]^^ ^ 1 +2k. The shear 
describes the area-preserving distortion of amplitude given by |y| and of direction given by 
its phase, see also Section WA\ and Fig. [151 In general the matrix Si also contains an anti- 
symmetric part associated with a rotation of the image but this term vanishes at linear order 
as can be seen from Eq. (|2.4l) . From Eq. (l2.6l) the convergence K and the shear components 
yi, y2, can be written at linear order in terms of the shear tensor as: 

^11+^22 ^11-^22 

K= y=yi + /y2 with yi= y2 = ^i2. (2.7) 

On the other hand, the gravitational potential, 4>, is related to the fluctuations of the density 
contrast, 6, by Poisson's equation: 

V^^^^-a^Hl{\+z)h with 5(x) = ^^4^, (2.8) 
1 P 

where p is the mean density of the universe. Note that since the convergence K and the 
shear components y,- can be expressed in terms of the scalar lensing potential (|) they are 
not independent. For instance, one can check from the first Eq. (l2.5l) and Eq. (12.71) that we 
have K 1 = yij +y2,2 (il49|) . T his a llows one to derive consistency relations satisfied by 



weak lensing distortions (e.g. (|242f) ) and deviations from these relations in the observed 



shear fields can be used to estimate the observational noise or systematics. Of course such 



relations also imply interrelations between correlation functions, see (12441) and Section [3] 
below. For a rigorous derivation of Eqs. (l2.4l) - (I2.7I) one needs to compute the paths of light 
rays (null geodesies) through the perturbed metric of spacetime using General Relativity 



152h . An altern ative approach is to follow the distortion of the cross-section of an infinites- 



imal light beam (|225l : l32l :l252l:ll7h. Both methods give back the results (l24l)-(l277]) obtained 
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in a heuristic manner above. 



2.2 Convergence, shear and aperture mass 



Thanks to the radial integration over % in (12.41) gradients of the gravitational potential along 
the radial direct ion give a negligible contribution as compared with transverse fluctuations 



(|146l : 11371 : 11701) since positive and negative fluctuations cancel along the line of sight. In 
other words , the radial integration selects Fourier radial modes of order \k\\ \ r^H /c (inverse 
of cosmological distances over which the effective lensing weight w{x) varies, see eg. (12. 101 ) 
below) whereas transverse modes are of order \k±\ ~ l/©0.5 ^ |^||| where 0.5 ^ 1 is the 
typical angular scale (a few arcmin) probed by the weak-lensing observable. Therefore, 
within this small-angle approximation the 2D Laplacian (12.71) associated with K can be 
expressed in terms of the 3D Laplacian (12.81) at each point along the line of sight. This 
yields for the convergence along a given line of sight up to Zs'- 

<zs) ^ Tdz w(x,%.)5(x) with w(x,%.) = 3^"^-^o'^(^)^(^---^) (i (2.9) 

Thus the convergence, k, can be expressed very simply as a function of the density field; 
it is merely an average of the local density contrast along the line of sight. Therefore, 
weak lensing observations allow one to measure the projected density field K on the sky 
(note that by looking at sources located at different redshifts one may also probe the radial 
direction). In practice the sources have a broad redshift distribution which needs to be taken 
into account. Thus, the quantity of interest is actually: 

K= / dz,n(z,)K(z,)= / dxwixMx) with h'(x) = / dz,n(z,) w(%,%,), (2-10) 

Jo Jo Jz 

where n{zs) is the mean redshift distribution of the sources (e.g. galaxies) normalized to 
unity and Zmux is the depth of the survey. Eg. (12. 101 ) neglects the discrete effects due to the 
finite number of galaxies, which can be obtained by taking into account the discrete nature 
of the distribution n(zs)- This gives corrections of order l/N to higher-order moments of 
weak-lensing observables, where A'^ is the number of galaxies within the field of interest. In 
practice is much larger than unity (for a circular window of radius 1 arcmin we expect 
> 100 for the SNAP mission) therefore it is usually sufficient to work with Eg. (12. 101) . 

In order to measure weak-lensing observables such as K or the shear y one measures for 
instance the brightness or the shape of galaxies located around a given direction on the 
sky. Therefore, one is led to consider weak-lensing guantities smoothed over a non-zero 
angular radius 0^ around the direction 0. More generally, one can define any smoothed 
weak-lensing guantity X(0) from its angular filter L'x(A0) by: 

X{Q) = I d0'[/x(0'-0)K(0O = I dxw I d0'C/x(0'-0)6(%,©0O, (2.11) 
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where 0' is the angular vector in the plane perpendicular to the line of sight (we restrict 
ourselves to small angular windows) and ©0' is the two-dimensional vector of transverse 
coordinates. Thus, it is customary to define the smoothed convergence by a top-hat L^k of 
angular radius 0^ but this quantity is not very convenient for practical purposes since it is 
easier to measure the ellipticity of galaxies (related to the shear y) than their magnification 
(related to k). This leads one to consider compensated filters Um^-p with polar symmetry 
which define the "aperture-mass" Map, that is with / dQUM^J^)— 0. Then Map can be ex- 
pressed in terms of the tangential component of the shear (|234|) so that it is not necessary 
to build a full convergence map from observations: 



Map(0) = I d0'f/M,p(|0'-0|)K(0') = I dd'QM,,{\Q'-m{Q') (2.12) 



where we introduce ( 234 ): 



Qm^p (0) = -f/M.p (0) + dQ' 0' Um^, (0') . 



(2.13) 



Besides, the aperture-mass provides a useful separation between E and B modes, as dis- 
cussed below in Section [33l 



For analytical and data analysis purposes it is often useful to work in Fourier space. Thus, 
we write for the 3D matter density contrast 6(x) and the 2D lensing potential (|)(0): 

5(x) = / ^^"''"S(k) and m = J j^e-''-'m, (2.14) 

where we use a flat-sky approximation for 2D fields. This is sufficient for most weak lensing 
purposes where we consider angular scales of the order of 1 — 10 arcmin, but we shall 
describe in Section l43] the more general expansion over spherical harmonics. From Eq. (|2.5l) 
and Pois son's equation (|2.8I) we obtain: 

«?) = -2/dxw(x)/^»-'"^p45,8(*||.4^;x), (2.15) 

where ^ is the component parallel to the line of sight of the 3D wavenumber k = {k^\,k±), 

with k± = i/T), and 5(k; %) is the matter density contrast in Fourier space at redshift z{x) ■ 
The weight w{x) along the line of sight was defined in Eqs. (|2.9l) - (I2.10I) . Then, from Eq. (l2.7l) 
we obtain for the convergence k: 

^i^) = -\i^l + ^'y)m - / /^^"''"'S [h^^y^) ■ (2.16) 

In the last expression we used as for Eq. (|2.9l ) Limber's approximation ~ k^ as the inte- 
gration along the line of sight associated with the projection on the sky suppresses radial 
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modes as compared with transverse wavenumbers (i.e. \k\\ \ <^ k^). In a similar fashion, we 
obtain from Eq. (|2.7l) for the complex shear y. 

m = -:^{(^. + i(^y)^m = ' /^^2 ^' ^^) = (2.17) 

X y 

where a is the polar angle of the wavenumber i = {(,x,(-y)- This expression clearly shows 
that the complex shear y is a spin-2 field: it transforms as y — >^ yg-^'^^ under a rotation of 
transverse coordinates axis of angle \|/. This comes from the fact that an ellipse transforms 
into itself through a rotation of 180 degrees and so does the shear which measures the 
area-preserving distortion, see Fig. [15] 

For smoothed weak-lensing observables X as defined in Eg. (12. Ill) we obtain: 

Xil) =Wx{-IQsW) with W^x(^e,) = y"d0/-^f/x(0), (2.18) 

where we introduced the Fourier transform Wx of the real- space filter Ux of angular scale 
0^. This gives for the convergence and the shear smoothed with a top-hat of angular radius 
0.: 



Wk(/0.O = WV(?0.) = Wk(£0.) (2.19) 



where J\ is the Bessel function of the first kind of order 1. In real space this gives back 
(with0= |0|): 

^.,e).^. u,,l).y-^,^, (2.20, 

where is the Heaviside function and (3 is the polar angle of the angular vector 0. Note 
that Eg. (12. 201) clearly shows that the smoothed convergence is an average of the density 
contrast over the cone of angular radius 0? whereas the smoothed shear can be written as 
an average of the density contrast outside of this cone. 

One drawback of the shear components is that they are even guantities (their sign can be 
changed through a rotation of axis, see Eg. (12. 171) ). hence their third-order moment vanishes 
by symmetry and one must measure the fourth-order moment (yf) (i.e. the kurtosis) in 
order to probe the deviations from Gaussianity. Therefore it is more convenient to use the 
aperture-mass defined in Eg. (12. 121) which can be derived from the shear but is not even, so 
that deviations from Gaussianity can be detected thro ugh the third-order moment (M^p). A 
simple example is provided by the pair of filters (234"): 

and: 



2474(^0.) 

WM,p(a,)= ^^Q^^2 • (2.22) 
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2.3 Approximations 



The derivation of Eq. (|2.4l) does not assume that the density fluctuations 5 are small but it 
assumes that deflection angles 50 are small so that the relative deflection of neighboring 
light rays can be computed from the gravitational potential gradients along the unperturbed 
trajectory (Born approximation). This may not be a good approximation for individual 
light beams, but in cosmological weak-lensing studies considered in this review one is 
only interested in the statistical properties of the gravitational lensing distortions. Since 
the statistical properties of the tidal field 4>,j are, to an excellent approximation, identical 
along the perturbed and unperturbed paths, the use of Eg. (12.41 ) is well-j ustifi ed to compute 



statistical quantities such as the correlation functions of the shear field (|146l : l22|) 



Apart from the higher-order corrections to the Born approximation discussed above (multi- 
ple lens couplings), other higher-order terms are produced by the observational procedure. 
Indeed, in Eg. (12. 101) we neglected the fluctuations of the galaxy distribution n{zs) which 
can be coupled to the matter density fluctuations along the line of sight. This source-lens 
correlation effect is more important as the overlapping area between the distributions of 
sources and lenses increases. On the other hand, source density fluctuations themselves can 
lead to spurious small-scale power (as the average distance to the sources can vary with the 
direction on the sky). Using analytical methods Ref. (l^ found that both these effects are 
negligible for the skewness and kurtosis of the convergence provided the source redshift 
dispersion is less than about 0.15. These source clustering effects were further discussed 



in (|93|) who found that numerical simulations agree well with semi-analytical estimates 



and that the amplitude of such effects strongly depends on the redshift distribution of the 



sources. A recent study of the source-lens clustering (1741) . using numerical simulations cou- 
pled to realistic semi-analytical models for the distribution of galaxies, finds that this effect 
can bias the estimation of Og by 2% — 5%. Therefore, accurate photometric redshifts will 
be needed for future missions such as SNAP or LSST to handle this effect. 
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3 Statistics of 2D Cosmic Shear 



For statistical analysis of cosmic shear, it is most common to use 2-point quantities, i.e. 
those which are quadratic in the shear, and calculated either in real or harmonic space. For 
this Section, we will restrict the discussion to 2D fields, where we consider the statistics 
of the shear pattern on the sky only, and not in 3D. The shear field will be treated as a 3D 
field in Section HI Examples of real-space 2-point statistics are the average shear variance 
and various shear correlation functions. In general there are advantages for cosmological 
parameter estimation in using harmonic- space statistics, as their correlation properties are 
more convenient, but for surveys with complicated geometry, such as happens with removal 
of bright stars and artifacts, there can be practical advantages to using real-space measures, 
as they can be easier to estimate. All the 2-point statistics can be related to the underlying 
3D matter power spectrum via the (2D) convergence power spectrum P^i^), and inspection 
of the relationship between the two point statistic and /'k(^) can be instructive, as it shows 
which wavenumbers are picked out by each statistic. In general, a narrow window in £ 
space may be desirable if the power spectrum is to be estimated. 



3. 1 Convergence and shear power spectra 

We define the power spectra P{k) of the 3D matter density contrast and /k(^) of the 2D 
convergence as: 

(5(ki)6(k2)) = (27r)36o(ki +V.2)P{ki) (3.1) 

and: 

WiWi)) = {2nfWi+h)PK{h). (3.2) 

The Dirac functions express statistical homogeneity whereas statistical isotropy implies 
that /'(k) and /k(^) only depend on k = |k| and £ = In Eq. (|3.2| ) we used a flat-sky 
approximation which is sufficient for most weak-lensing purposes. We shall discuss in 
Section |4] the expansion over spherical harmonics (instead of plane waves as in Eq. (|3.2|) ) 
which is necessary for instance for full-sky studies. Then, from Eqs. (l2.16l) - (I2.17I) we obtain: 

P^{^)^P^{l)=^-tp^{l) (3.3) 

and: 

Thus this expression gives the 2D convergence power spectrum in terms of the 3D matter 
power spectrum P{k;x) integrated along the line of sight, using Limber's approximation. 
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3.2 2-point statistics in real space 



As an example of a real-space 2-point statistic, consider the shear variance, defined as 
the variance of the average shear y evaluated in circular patches of varying radius G.,. The 
averaging is a convolution, so the power is multiplied (see Eqs. (l2.18l) , (|2.19l )): 

where Jn is a Bessel function of order n. 

The shear correlation functions can either be defined with reference to the coordinate axes, 

^,,(0)^(y,(eOy,(e'+0)) (3.6) 

where /, j = 1,2 and the averaging is done over pairs of galaxies separated by angle = |0| . 
By parity ^12 = 0, and by isotropy and ^22 are functions only of |0|. The correlation 
function of the complex shear is 



(Yy*)e = /^W/-^^ 0.1) 



I 



f/if 



(271)2 

= / ^iP^{e)Jom- 

Alternatively, the shears may be referred to axes oriented tangentially (t) and at 45 degrees 
to the radius (x), defined with respect to each pair of galaxies used in the averaging. The 
rotations y^ y = ye^^'^, where \|/ is the position angle of the pair, give tangential and cross 
components of the rotated shear as Y = — jt — ijx , where the components have correlation 
functions ^tt and x respectively. It is common to define a pair of correlations 

^±(0)=^„±^xx, (3.8) 



which can be related to the convergence power spectrum by (see (|146|) ) 



r°° df 

^_(0)=/ ip^{i)j,{m). (3.9) 

Finally, let us consider the class of statistics referred to as aperture masses associated with 
compensated filters, y^hi ch. we defined in equation (12.121) . This allows Map to be related to 



the tangential shear (|234r ) as in equation (I2.12|) . Several forms of Um^^ have been suggested. 
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Fig. 2. Kernel functions for the two-point statistics discussed in this Section, z = iOs- The thin solid 
line peaking at z = is corresponds with the shear variance, the thick solid line is the aperture mass, 
with filter given in equation (12.21b . the short dashed line is the kernel for ^+ and the long dashed to 

which trade locality in real space with locality in i space. Ref. considers the filter 

f/^fap of equation (12.211) which cuts off at some scale, 0^. From equation (12.221) this gives a 
two-point statistic 

Other forms have been suggested ^Ed), which are broader in real space, but pick up a nar- 
rower range of i power for a given 0. As we have seen, all of these two-point statistics can 
be written as integrals over i of the convergence power spectrum / k(^) multiplied by some 
kernel function, since weak-lensing distortions can be expressed in terms of the lensing 
potential (|), see Eqs. (|2.5l) and (12.161) . 

If one wants to estimate the matter power spectrum, then there are some advantages in hav- 
ing a narrow kernel function, but the uncertainty principle then demands that the filtering 
is broad on the sky. This can lead to practical difficulties in dealing with holes, edges etc. 
Filter functions for the 2-point statistics mentioned here are shown in Fig. [21 

3.3 E/B decomposition 

Weak gravitational lensing does not produce the full range of locally linear distortions 
possible. These are characterised by translation, rotation, dilation and shear, with six free 
parameters. Translation is not readily observable, but weak lensing is specified by three 
parameters rather than the four remaining degrees of freedom permitted by local affine 
transformations. This restriction is manifested in a number of ways: for example, the trans- 
formation of angles involves a 2 x 2 matrix which is symmetric, so not completely general. 
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Fig. 3. Illustrative E and B modes: the E modes show what is expected around overdensities (left) 
and underdensities (right). The B mode patterns should not be seen (from van Waerbeke & Mellier 
2003). 

see equation (12.61) . Alternatively, a general spin-weight 2 field can be written in terms of 
second derivatives of a complex potential, whereas the lensing potential is real. As noticed 
below equation (12.81) and in equation (13.251) . this also implies that there are many other 
consistency relations which have to hold if lensing is responsible for the observed shear 
field. In practice the observed ellipticity field may not satisfy the expected relations, if it is 
contaminated by distortions not associated with weak lensing. The most obvious of these 
is optical distortions of the telescope system, but could also involve physical effects such 
as intrinsic alignment of galaxy ellipticities, which we will consider in Section|4l 

A convenient way to characterise the distortions is via E/B decomposition, where the shear 
field is described in terms of an 'E-mode', which is allowed by weak lensing, and a 'B- 
mode', which is not. These terms are borrowed from similar decompositions in polarisation 
fields. In fact weak lensing can generate B-modes, but they are expected to be very small 



(|240|) . so the existence of a significant B-mode in the observed shear pattern is indicative 
of some non-lensing contamination. Illustrative examples of E- and B-modes are shown in 
Fig|3](from('299)). The easiest way to introduce a B-mode mathematically is to make the 
lensing potential complex: 

(^ = (^£ + %. (3.11) 

There are various ways to determine whether a B-mode is present. A neat way is to gen- 
eralise the aperture mass to a complex M = Map -|- iM^, where the real part picks up th e E 
modes, and the imaginary part the B modes. Alternatively, the ^± can be used ([60; 241 ): 



poo 

Pyc±{i)=Ti je0[7o(^e)^+(0)±74(^e)^-(0)] (3.12) 

Jo 

where the ± power spectra refer to E and B mode powers. In principle this requires the 
correlation functions to be known over all scales from to °o. Variants of this (|60l) allow 
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the E/B-mode correlation functions to be written in terms of integrals of ^± over a finite 
range: 



^£(0) = i[^_(0)+^V(e)] 
Ue) = -^[^-(e)-^V(0 



(3.13) 



where 



This avoids the need to know the correlation functions on large scales, but needs the ob- 
served correlation functions to be extrapolated to s mall scales; this was one of the ap- 
proaches taken in the analysis of the CFHTLS data (Ill7|) . Difficulties with estimating the 
correlation functions on small scales ha ve le d others to prefer to extrapolat e to l arge scales, 
such as in the analysis of the GEMS (jlOlb and William Herschel data (jlSSb . Note that 
without full sky coverage, the decomposition into E and B modes is ambiguous, although 
for scales much smaller than the survey it is not an issue. 



3.4 Estimators and their covariance 



The most common estimate of the cosmic shear comes from measuring the ellipticities of 
individual galaxies. We will consider the practicalities in Section [6l For weak gravitational 
lensing, these estimates are very noisy, since the galaxies as a population have intrinsic 
ellipticities es with a dispersion of about 0.4, whereas the typical cosmic shear is around 
Y~ 0.01. Therefore, one needs a large number of galaxies to decrease the noise ~ es/VN 
associated with these intrinsic ellipticities (hence one needs to ob serv e the more numerous 



faint galaxies). The observed ellipticity is related to the shear by ( 233b 



e — 2g-\-g^e* 

where ^ = y/ ( 1 — k) is the reduced shear. Here we defined the ellipticity e such that for an 
elliptical image of axis ratio r < 1 we have: 



l-r2 



l+r 



,2' 



(3.16) 



Other definitions are also used in the literature such as = (1 — r)/(l + r), see (|244|) . To 
linear order in y or K, we obtain from Eq. (|3.15l) : 

e~e5 + 2y, (3.17) 

with a small correction term when averaged, e is therefore dominated by the intrinsic ellip- 
ticity, and many source galaxies are needed to get a robust measurement of cosmic shear. 
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This results in estimators of averaged quantities, such as the average shear in an aperture, 
or a weighted average in the case of Map. Any analysis of these quantities needs to take 
account of their noise properties, and more generally in their covariance properties. We 
will look only at a couple of examples here; a more d etaile d discussion of covariance of 
estimators, including non-linear cumulants, appears in (|202|) . 



3.4.1 Linear estimators 

Perhaps the simplest average statistic to use is the average (of A'^) galaxy ellipticities in a 
2D aperture on the sky: 

Y4|f. (3.18) 
The covariance of two of these estimators ya and is 

1 N M 

(TocYp) = i;((^« + 2y,)(4, + 2YP), (3.19) 

where the apertures have and M galaxies respectively. If we assume (almost certainly 
incorrectly; see Section 14.71) that the source ellipticities are uncorrelated with each other, 
and with the shear, then for distinct apertures the estimator is an unbiased estimator of the 
shear correlation function averaged over the pair separations. If the apertures overlap, then 
this is not the case. For example, in the shear variance, the apertures are the same, and 

(It'|) = ^II(^5,,+y,^), (3.20) 

which is dominated by the presence of the intrinsic ellipticity variance, = (1^5^) ~ 
0.3^-0.42. The average shear therefore has a variance of ol/ AN. If we use the (quadratic) 
shear variance itself as a statistic, then it is estimated by omitting the diagonal terms: 

1 ^ 

\f\ = ^ ^Y.T.e^e). (3.21) 

" ' 4A^(A^-1) f^^j^. J ^ ^ 

For aperture masses (equation (12.121) . the intrinsic ellipticity distribution leads to a shot 
noise term from the finite number of galaxies. Again we simplify the discussion here by 
neglecting co rrela tions of source ellipticities. The shot noise can be calculated by the stan- 
dard method (Si of dividing the integration solid angle into cells i of size A G, containing 
Hi = or 1 galaxy: 

Map - £A20,n,2(|e,|) {es,/2 + y,)t. (3.22) 

Squaring and taking the ensemble average, noting that {esijesj^t) = <5eS/j/2, nf = ni, and 
rewriting as a continuous integral gives 

{M%)sN^^ f d^QQWQ\). (3.23) 
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Shot noise terms for other statistics are calculated in similar fashion. In addition to the 
covariance from shot noise, there can be signal covariance, for example from samples of 
different depths in the same area of sky; b oth s amples are affected by the lensing by the 
common low-redshift foreground structure (|203|) . 



3.4.2 Quadratic estimators 

We have already seen how to estimate in an unbiased way the shear variance. The shear 
correlation functions can similarly be estimated: 

^,(9) = ^'^"'''"^<''-"^'^^'"^^"^ (3.24) 

^Y.ijWiWj 

where the w,- are arbitrary weights, and the sum extends over all pairs of source galaxies 
with separations close to 0. Only in the absence of intrinsic correlations, {euep ± e/x = 
OgS/j + 4^±(|0/ — 0j|), are these estimators unbiased. The variance of the shear (|yP) and 
of the aperture-mass (M^p) can also be obtained from the shear correlation functions (as 
may be seen for instance from Eg. (13. 101) ). This avoids the need to place circular apertures 
on the sky which is hampered by the gaps and holes encountered in actual weak lensing 
surveys. 

As with any quadratic quantity, the covariance of these estimators depends on the 4-point 
function of the source ellipticities and the shear. These expressions can be evaluated if 
the shear field is assumed to be Gaussian, but the expressions for this (and the squared 
apert ure mass covariance) are too cumbersome to be given here, so the reader is directed to 
( 24lb . At small angular scales (below ~ 10'), which are sensitive to the non-linear regime of 



gravitational clustering, the fields can no longer be approximated as Gaussian and one must 
use numerica l sim ulations to calibrate the non-Gaussian contributions to the covariane, as 



described in (|258r) . 



In harmonic space, the convergence power spectrum may be estimated from either ^+ or 
^_ (or both), using Eq. (13.91) . From the orthonormality of the Bessel functions, 

P^{e)= / rf00^±(0)7o,4(^0) (3.25) 
Jo 

where the 0, 4 correspond to the +/— cases. In practice, ^±(0) is not known for all 0, and 
the integral is truncated on both small and large scales. This can lead to inaccuracies in 
the estimation of /'k(^) (see (241)). An alternative method is to parametrise /k(^) in band- 
powers, a nd to use parameter estimation techniques to estimate it from the shear correlation 



functions (1231:1381 
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3.4.3 2-point statistics measurement 



Since the first measurements of weak lensing by large scale structures (l295l:[6l: [309l : [l54l) . 



all ideas discussed above have been put in practice on real data. Figure |4] shows the mea- 
surement of shear top-hat variance as function of scale. Some groups performed a E/B sep- 
aration, which lead to a more accurate measurement of residual systematics. Table H] shows 
all measurements of t he mass power spectrum as to date. It is interesting to note, except 



for COMBO- 17 (|101|) . none of the measurements are using a source redshift distribution 



obtained from the data, the y all use the Hubble Deep Fields with different prescriptions 



regarding galaxy weighing. (|300t) have shown that the Hubble Deep Field photometric red- 
shift distribution can lead to a source mean redshift error of ~ 10% due to cosmic variance. 
The relative tension between different measures of Og in Table |4] comes in part from this 
problem, and also from an uncertainty regardin g how to treat the residual systematics, the 



ngnov 
■ 1l33 



5-mode, in the cosmic shear signal. Recently, (11331) have released the largest photomet- 
ric redshift catalogue, obtained from the CFHTLS-DEEP data. The most recent 2-point 
statitics analysis involves the combination of this photometric redshift sample with the 
largest weak lensing surveys described in Table IH In this analysis the relative tension is 
gone and points towards a value of the power spectrum normalisation Og = 0.75 ±0.05 if 
one considers all possible issues with uncertainties in the source redshift distribution (20). 
From Table |4] and Figure |4] one can say that the amplitude of weak lensing by large scale 
structure was successfully measured and the main uncertainty remains the source redshift 
distribution. The systematics due to the Point Spread Function correction (see Section (6]) 
seems to be much better understood than a few years ago, and many promising techniques 
have been proposed to solve it. The proper calibration of the redshift distribution is likely 
to remain the major limitation for the use of weak lensing in precision cosmology. 
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Fig. 4. Compilation of most of the shear measurements Usted in Tabled The vertical axis is the shear 
top-hat variance multiplied by the angular scale in arcminutes. The horizontal axis is the radius of 
the smoothing window in arcminutes. The positioning along the y-axis is only approximate given 
that the different surveys have a slightly different source redshift distribution. The RCS result (mean 
source redshift of 0.6) was rescaled to a mean source redshift of one. 
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Table 4 

Reported constraints on the power spectrum normalization "ag" for = 0.3 for a flat Universe, obtained from a given "statistic" (from and 
extended). "CosVar" tells us whether or not the cosmic variance has been included, "E/B" tells us whether or not a mode decomposition has been 
used in the likelihood analysis. z.s and F are the priors used for the different surveys identified with "ID". 



ID 


08 (n = 0.3) 


Statistic 


Field 




CosVar 


E/B 




r 


Maoh et al. 01 


1.03 ±0.05 


if) 


VLT+CTIO+WHT+CFHT 


- 


no 


no 


- 


0.21 


Van Waerbeke et al. 01 


0.88±0.11 


{f),m{Mi) 


CFHT 8 sq.deg. 


1=24.5 


no 


no (yes) 


1.1 


0.21 


Rhodes et al. 01 






HST 0.05 sq.deg. 


1=26 


yes 


no 


0.9-1.1 


0.25 


Hoekstra et al. 02 


0.81 ±0.08 


if) 


CFHT+CTIO 24 sq.deg. 


R=24 


yes 


no 


0.55 


0.21 


Bacon et al. 03 


0.97 ±0.13 


k{r) 


Keck+WHT 1.6 sq.deg. 


R=25 


yes 


no 


0.7-0.9 


0.21 


Refregier et al. 02 


0.94 ±0.17 


if) 


HST 0.36 sq.deg. 


1=23.5 


yes 


no 


0.8-1.0 


0.21 


Van Waerbeke et al. 02 


0.94 ±0.12 


{Ml,) 


CFHT 12 sq.deg. 


1=24.5 


yes 


yes 


0.78-1.08 


0.1-0.4 


Hoekstra et al. 02 


n Q 1+0 05 


{f),m{Mi,) 


CFHT+CTIO 53 sq.deg. 


R=24 


yes 


yes 


0.54-0.66 


0.05-0.5 


Brown et al. 03 


0.74 ±0.09 




C0MB017 1.25 sq.deg. 


R=25.5 


yes 


no (yes) 


0.8-0.9 


- 


Hamana et al. 03 


(2a)0.69l°i5 


{Ml,),k{r) 


Subaru 2.1 sq.deg. 


R=26 


yes 


yes 


0.8-1.4 


0.1-0.4 


Jarvis et al. 03 


(2a)0.7ll«:}^ 


in m {Ml,) 


CTIO 75 sq.deg. 


R=23 


yes 


yes 


0.66 


0.15-0.5 


Rhodes et al. 04 


1.02±0.16 


{'f),Ur) 


STIS 0.25 sq.deg. 


(/) = 24.8 


yes 


no 


1.0 ±0.1 




Hey mans et al. 05 


0.68±0.13 


if), m 


GEMS 0.3 sq.deg. 


(wJeoe) =25.6 


yes 


no (yes) 


~ 1 




Massey et al. 05 


1.02±0.15 


if), m 


WHT 4 sq.deg. 


R=25.8 


yes 


yes 


~0.8 




Van Waerbeke et al. 05 


0.83 ±0.07 


{'f),m 


CFHT 12 sq.deg. 


1=24.5 


yes 


no (yes) 


0.9 ±0.1 


0.1-0.3 


Heitterscheidt et al. 06 


0.8 ±0.1 


if), m 


GaBoDS 13 sq.deg. 


R=[21. 5,24.5] 


yes 


yes 


~0.78 


he [0.63,0.77] 


Semboloni et al. 06 


0.90±0.14 


{Ml,), m 


CFHTLS-DEEP 2.3 sq.deg. 


i=25.5 


yes 


yes 


~ 1 


T = ah 


Hoekstra et al. 06 


0.85 ±0.06 


{f),%{r),{Ml,) 


CFHTLS-WIDE 22 sq.deg. 


i=24.5 


yes 


yes 


0.8 ±0.1 


r = D.h 



3.5 Mass Reconstruction 



The problem of mass reconstruction is a central topic in weak lensing. Historically this is 
because the early measurements of weak gravitational lensin g we re obtained in clusters 



of galaxies, and this led to the very first maps of dark matter (|33l : 1731) . These maps were 
the very first demonstrations that we could see the dark side of the Universe without any 
assumption regarding the light-mass relation, which, of course, was a major breakthrough 
in our exploration of the Universe. Reconstructing mass maps is also the only way to per- 
form a complete comparison of the dark matter distribution to the Universe as seen in other 
wavelengths. For these reasons, mass reconstruction is also part of the shear measurement 
process. A recent example of the power of mass reconstruction is shown by the bullet clus- 



ter (1471: 134|) . which clearly indicates the presence of dark matter at a location different from 
where most of the baryons are. This is a clear demonstration that, at least for the extreme 
cases where light and baryons do not trace the mass, weak lensing is the only method that 
can probe the matter distribution. 

A mass map is a convergence, K map (projected mass), which can be reconstructed from 
the shear field y, : 



K = i {dl + aj) ; Yi = i {dl - a2) ; 72 = a,-a,(^, (3.26) 

where ^ is the projected gravitational potential (146'), see Eqs. (l2.5l) - (|2.7I) . Assuming that 
the reduced shear and shear are equal to first approximation, i.e. gi ~ y, (which is true only 
in the weak lensing regime when |y| ^ 1 and K -C 1), (fl47) have shown that the Fourier 
transform of the smoothed convergence k{£) can be obtained from the Fourier transform of 
the smoothed shear map y{i) : 

This follows from Eq. (|2.17l) if smoothing is a mere convolution as in Eq. (|2.1 II ) which writes 
in Fourier space as a mere product, see Eq. (l2.18l) . This relation explicitly shows that mass 
reconstruction must be performed with a given smoothing window, otherwise the variance 



of the mass map becomes infinite (Il47t ). Indeed, the random galaxy intrinsic ellipticities est 
introduce a white noise which gives a large-£ divergence when we transform back to real 
space for the variance (k^)c. The Fourier transform method is a fast NlogN process, but the 
non-linear regions K ~ y, ~ 1 are not accurately reconstructed. A likelihood reconstruction 
method works in the intermediate and strong lensing regimes (P). A function of the 
reduced shear gi is minimised by finding the best gravitational potential ^ij calculated on a 
grid ij: 



"7 



(3.28) 
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Fig. 5. Left: simulated noise-free K map. Field-of-view is 49 square degrees in LCDM cosmology. 
Right: reconstructed K field with reaUstic noise level (give details). Van Waerbeke et al. (1999) have 
shown that two and three-point statistics can be accurately measured from such mass reconstruc- 
tions. 



The shot noise of the reconstructed mass map depends on the sm oothi ng window, the in- 
trinsic ellipticity of the galaxies and the number density of galaxies (|294r) . The tw o met hods 
outlined above provide an accurate description of the shot noise: it was shown (|293l) that 
two and three-points statistics can be measured accurately from reconstructed mass maps 
using these methods (see Fig ure | 5l). Im porta nt for cluster lensing, the non-linear version of 
([147) has been developed in (Il5l|y . and ( 233 ) have developed an alternative which also con- 
serves the statistical properties of the noise. The advantage of a reconstruction method that 
leaves intact the shot noise is that a stati stical analysis of the mass map is relatively straight- 
forward (e.g. the peak statistics in (S))._ Mass reconstruction has proven to be reasonably 
successful in blind cluster searches (1981 : Ssl; ItSi) a radically different approach in map 
making consists in reducing the noise in ord er to iden tify t he highest signal- to-no ise peaks. 
Such an approach has been developed by (12531 : l36l : |184|) . More recently proposed 
a wavelet approach, where the size of the smoothing kernel is optimized as a function of 
the local noise amplit ude. An application of this method on the COSMOS data for cluster 
detection is shown in ( 187 ). 



Mass maps are essential for some specific cosmological stu dies s uch as morphology anal- 
ysis like Minkowski functionals and Euler characteristics (228) and for global statistics 
(Probability distribution function of the convergence, e.g. (316). The reconstruction pro- 
cesses is non local, and this is why it is difficult to have a perfect control of the error 
propagation and systematics in the K maps. In particular, note that one can only reconstruct 
the convergence K up to a constant Kq, since Eg. (13. 271) is undetermined for i = 0. Indeed, 
a constant matter surface density in the lens plane does not create any shear (since it se- 
lects no prefered direction) but it leads to a non-zero constant convergence (which may 
only be eliminated if we observe a wide-enough field where the mean convergence should 
vanish, or use complementary information such as number counts which are affected by 
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the associated mag nification) . This is the well-known mass-sheet degeneracy. The aperture 



mass statistics Map (| 1481 : 12371) introduced in Eg. (12. 121) has been invented to enforce locality 



of the mass reconstruction, therefore it might provide an alternative to the inversion prob- 
lem, although it is not yet clear that it can achieve a signal-to-noise as good as top-hat or 
gaussian smoothing windows. 

The aperture mass statistic c an also be used to provide unbiased estimates of the power 
spectrum QS), galaxy biasing (12921 : |228|) . high order statistics (12371 : |296|) . and peak statistics 
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4 3D Weak Lensing 



4.1 What is 3D weak lensing? 

The way in which weak lensing surveys have been analysed to date has been to look for 
correlations of shapes of galaxies on the sky; this can be done even if there is no distance 
information available for individual sources. However, as we have seen, the interpretation 
of observed correlations depends on where the imaged galaxies are: the more distant they 
are, the greater the correlation of the images. One therefore needs to know the statistical 
distribution of the source galaxies, and ignorance of this can lead to relatively large errors 
in recovered parameters. In order to rectify this, most lensing surveys obtain multi-colour 
photometry of the sources, from which one can estimate their redshifts. These 'photometric 
redshifts' are not as accurate as spectroscopic redshifts, but the typical depth of survey 
required by lensing surveys makes spectroscopy an impractical option for large numbers 
of sources. 3D weak lensing uses the distance information of individual sources, rather 
than just the distribution of distances. If one has an estimate of the distance to each source 
galaxy, then one can utilise this information and investigate lensing in three dimensions. 
Essentially one has an estimate of the shear field at a number of discrete locations in 3D. 

There are several ways 3D information can be used: one is to reconstruct the 3D grav- 
itational potential or the overdensity field from 3D lensing data. We will look at this in 
Section 14.21 The second is to exploit the additional statistical power of 3D information, 
firstly by dividing the sources into a number of shells based on estimated redshifts. One 
then essentially performs a standard lensing analysis on each shell, but exploits the ex- 
tra information from cross-correlations between shells. This sort of analysis is commonly 
referred to as tomography, and we explore this in Section 14.31 and in Section 14.41 where 
one uses ratios of shears behind clusters of galaxies. Finally, one can perform a fully-3D 
analysis of the estimated shear field. Each approach has its merits. We cover 3D statistical 
analysis in Section |431 At the end of Section |477l we investigate how photometric redshifts 
can remove a potentially important physical systematic: the intrinsic alignment of galaxies, 
which could be wrongly interpreted as a shear signal. Finally, in Section [4^81 we consider a 
potentially very important systematic error arising from a correlation between cosmic shear 
and the intrinsic alignment of foreground galaxies, which could arise if the latter responds 
to the local tidal gravitational field which is partly responsible for the shear. 



4.2 3D potential and mass reconstruction 

As we have already seen, it is possible to reconstruct the surface density of a lens system by 
analysing the shear pattern of galaxies in the background. An interesting question is then 
whether the reconstruction can be done in three dimensions, when distance information 
is available for the sources. It is probably self-evident that mass distributions can be con- 
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strained by the shear pattern, but the more interesting possibility is that one may be able to 
determine the 3D mass density in an essentially non-parametric way from the shear data. 



The idea (|280l) is that the shear pattern is derivable from the lensing potential (|)(r), which 
is dependent on the gravitational potential ^'(r) through the integral equation 



*w 41 "''(7 -7) *(■■'' 



(4.1) 



where the integral is understood to be along a radial path (the Bom approximation), and 
a flat Universe is assumed in equation (14.11) . In this Section we work with spherical co- 
ordinates (we do not use the flat-sky approximation except where explicitely noticed) and 
we note r = rn the position on the line of sight in the direction n, at the comoving radial 
distance r (which we also noted % in Section [21). The gravitational potential is related to 
the density field via Poisson's equation (12.81) . There are two problems to solve here; one 
is to construct (|) from the lensing data, the second is to invert equation (14.11) . The second 
problem is straightforward: the solution is 



4>(r) 



c2 a_ 

2 dr 



r^-^ir) 



(4.2) 



From this and Poisson's equation = {3/2)HQQ.^d/a{t), we can reconstruct the mass 
overdensity field 



6(r) 



a{t)c^ 

3//o^ni 



r^-m 



(4.3) 



The construction of (|) is more tricky, as it is not directly observable, but must be estimated 
from the shear field. This reconstruction of the lensing potential suffers from a similar 
ambiguity to the mass-sheet degeneracy for simple lenses. To see how, we first note that 
the complex shear field y is the second derivative of the lensing potential (Eq. (|2.7l) ): 



y(r) 
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1/a^ 

2 \dx^ dy^ J dxdy 



(4.4) 



As a consequence, since the lensing potential is real, its estimate is ambiguous up to the 
addition of any field /(r) for which 



a2/(r) a2/(r) a2/(, 



0. 



3^2 9j2 dxdy 
Since (|) must be real, the general solution to this is 

/(r) =F{r)+ G{r)x + H{r)y + P{r) {x^ +y^) 



(4.5) 



(4.6) 



where F, G, H and P are arbitrary functions of r = |r|. Assuming these functions vary 
smoothly with r, only the last of these survives at a significant level to the mass density 
(since the 3D Laplacian V2 in Eq. (|4.3l) is dominated by the 2D Laplacian d^ + dy as for 
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Eq. (|2.9l )). and it corresponds to a sheet of overdensity 



4a{t)c^ d 



(4.7) 



There are a couple of ways to deal with this problem. For a reasonably large survey, one 
can assume that the potential and its derivatives are zero on average, at each r, or that the 
overdensity has average value zero. For further details, see (8). Note that the relationship 
between the overdensity field and the lensing potential is a linear one, so if one chooses a 
discrete binning of the quantities, one can use standard linear algebra methods to attempt 
an inversion, subject to some constraints such as minimising the expected recon struc tion 
errors. With prior knowledge of the signal properties, this is the Wiener filter. See (|l28h for 
further details of this approach. 



4.3 Tomography 



In the case where one has distance information for individual sources, it makes sense to em- 
ploy the information for statistical studies. A natural course of action is to divide the survey 
into slices at different distances, and perform a study of the shear pattern on each slice. In 
order to use the information effectively, it is necessary to look at cross -correlations of the 



shear fields in the slices, as well as correlations within each slice (Ill9|) . This procedure is 



usually referred to as tomography, although the term does not seem entirely appropriate. 

We start by considering the average shear in a shell, which is characterised by a probability 
distribution for the source redshifts z = z{r), p{z) ■ We also take the opportunity to introduce 
the more complicated edth derivative ( 9 ) on the curved sky, which is required if the survey 
does not subtend a small angle on the sky We shall not use make use of this much in 
this review, but it is included for completeness. The shear is the second edth derivative of 
the lensing potential , 

Y(r) = ^5 g(^(r) ^ l(a, + /a,)2(^(r) (4.8) 

where the last equality holds in the flat-sky limit. If we average the shear in a shell, giving 
equal weight to each galaxy, then the average shear can be written in terms of an effective 
lensing potential 

poo 

(^eff(e) = / dzp{z)m (4.9) 

where the integral is at fixed 0, and p{z) is zero outside the slice (we ignore errors in 
distance estimates such as photometric redshifts; these could be incorporated with a suitable 
modification to p{z))- In terms of the gravitational potential, the effective lensing potential 
is 

c^eff(e) = -/ dr^{r)g{r) (4.10) 
Jo 
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where reversal of the order of integration gives the lensing efficiency to be 



dz p{z) 



1 1 



(4.11) 



where z! = z'{r') and we assume flat space. If we perform a spherical harmonic transform 
of the effective potentials for slices i and j, then the cross power spectrum can be related to 
the power spectrum of the gravitational potential P(^{k) via a version of Limber's equation: 



(4.12) 



where 



(4.13) 



is the cross power spectrum of the lensing potentials. The last argument in allows for 
evolution of the power spectrum with time, or equivalently distance. The power spectra of 
the convergence and shear are related to ctf ■ by ( 120 ) 



CflJ-^^C% (4.14) 

_l (£ + 2)! 

4(£-2)! 

The sensitivity of the cross power spectra to cosmological parameters is through various 
effects, as in 2D lensing: the shape of the linear gravitational potential power spectrum is 
dependent on some parameters, as is its nonlinear evolution; in addition the z{r) relation 
probes cosmology. The reader is referred to standard cosmological texts for more details of 
the dependence of the distance-redshift relation on cosmological parameters. 



Ref. (|ll9l) illustrates the power and limitation of tomography, with two shells (Fig. (6]). As 
expected, the deeper shell (2) has a larger lensing power spectrum than the nearby shell (1), 
but it is no surprise to find that the power spectra from shells are correlated, since the light 
from both passes through some common material. Thus one does gain from tomography, 
but, depending on what one wants to measure, the gains may or may not be very much. 
For example, tomography adds rather little to the accuracy of the amplitude of the power 
spectrum, but far more to studies of dark energy properties. One also needs to worry about 
systematic effects, as leakage of galaxies from one shell to another, throu gh noisy o r biased 
photometric redshifts, can degrade the accuracy of parameter estimation ( 132l : ll75 ). Figure 
|7] shows the first tentative of tomographic measurement using photometric redshift on the 
CFHTLS deep data. 
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Fig. 7. Tomographic shear measurement on the CFHTLS-deep (from Semboloni et al 2006). Yellow 
lines show the it la contours around high and low redshift lensed galaxies. The filled triangles and 
circles show the E-mode and the empty symbols the B-mode. 

4.4 The Shear Ratio test 



The shear contributed by the general large-scale structure is typically about 1%, but the 
shear behind a cluster of galaxies can far exceed this. As always, the shear of a background 
source is dependent on its redshift, and on cosmology, but also on the mass distribution in 
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the cluster. This can be difficuk to model, so it is attractive to consider methods which are 
decoupled fror n the det ails o f the mass distribution of the cluster. Various methods have 
been proposed (fl39l: l3d : 13151) . The method currently receiving the most attention is simply 



to take ratios of average tangential shear in different redshift slices for sources behind the 
cluster. 

The amplitude of the induced tangential shear is dependent on the source redshift Zs, and 



on cosmology via the angular diameter distance-redshift relation 'D [x{zs)] by (|28l|) : 



JAZ) JAz ) ^^^^^^^^ , C4.i3; 

where jt ^oo is the shear which a galaxy at infinite distance would experience, and which 
characterises the strength of the distortions induced by the cluster, at redshift zi- Evidently, 
we can neatly eliminate the cluster details by taking ratios of tangential shears, for pairs of 
shells in source redshift: 



^ _ lt,i ^ 'D[x{zj)]'D[x{z,)-x{zi)] 
" %j 'D[x{zi)]-D[x{zj)-X{zi)]- ^ • ^ 



In reality, the light from the more distant shell passes through an extra pathlength of clumpy 
matter, so suffers an additional source of shear. This can be treated as a noise term (281). 
This approach is attractive in that it probes cosmology through the distance-redshift relation 
alone, being (at least to good approximation) independent of the growth rate of the fluc- 
tuations. Its dependence on cosmological parameters is therefore rather simpler, as many 
parameters (such as the amplitude of matter fluctuations) do not affect the ratio except 
through minor side-effects. More significantly, it can be used in conjunction with lensing 
methods which probe both the distance-redshift relation and the growth-rate of structure. 
Such a dual approach can in principle distinguish between quintessence-type dark energy 
models and modifications of Einstein gravity. This possibility arises because the effect on 
global properties (e.g. z(%)) is different from the effect on perturbed quantities (e.g. the 
growth rate of the power spectrum) in the two cases. The method has a signal-to-noise 
which is limited by the finite number of clusters which are massive enough to have mea- 
surable tangential shear. In an all-sky survey, the bulk of the signal would come from the 
10^ — 10^ clusters above a mass limit of lO^^M©. 



4.5 Full 3D analysis of the shear field 



An alternative approach to take is to recognise that, with photometric redshift estimates 
for individual sources, the data one is working with is a very noisy 3D shear field, which 
is sampled at a number of discrete locations, and for whom the locations are somewhat 
imprecisely known. It makes some sense, therefore, to deal with the data one has, and to 
compare the statistics of the discrete 3D field with theoretical predictions. This was the 
approach of (|97|; |4l|; 1981) . It should yield smaller statistical errors than tomography, as it 
avoids the binning process which loses information. 
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In common with many other methods, one has to make a decision whether to analyse the 
data in configuration space or in the spectral domain. The former, usually studied via cor- 
relation functions, is advantageous for complex survey geometries, where the convolution 
with a complex window function implicit in spectral methods is avoided. However, the 
more readily computed correlation properties of a spectral analysis are a definite advantage 
for Bayesian parameter estimation, and we follow that approach here. 

The natural expansion of a 3D scalar field which is derived from a potential is in terms 
of products of spherical harmonics and spherical Bessel functions, j£{kr)Yp{B,(p), using 
the spherical coordinates (r, 0,(p). Such products, characterised by 3 spectral parameters 
{k,i,m), are eigenfunctions of the Laplace operator, thus making it very easy to relate the 
expansion coefficients of the density field to that of the potential (essentially via —k^ from 
the operator). Similarly, the 3D expansion of the lensing potential, 

<^e^{k) = \j\f d\^{r)kk{kr)Y^{Q,i^), (4.17) 

where the prefactor and the factor of k are introduced for convenience. The expansion of the 
complex shear field is most naturally made in terms of spin-weight 2 spherical harmonics 
2^™ and spherical Bessel functions, since 7=^9 5(|), and 3 57"' 2^7"' 

"^W = V / dkye^kMkr) 2ir(e,(p). (4.18) 

' £m 

The choice of the expansion becomes clear when we see that the coefficients of the shear 
field are related very simply to those of the lensing potential: 




(£ + 2)^ 

JlZ^^'^Uk). (4.19) 

The relation of the <^im{k) coefficients to the expansion of the density field is readily com- 
puted, but more complicated as the lensing potential is a weighted integral of the gravi- 
tational potential. The details will not be given here, but relevant effects such as photo- 
metric redshift errors, nonlinear evolution of the power spectrum, and the discreteness of 
the sampling are easily included. The reader is referred to the original papers for details 
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In this way the correlation properties of the Jimik) coefficients can be related to an integral 
over the power spectrum, involving the z{r) relation, so cosmological parameters can be 
estimated via standard Bayesian methods from the coefficients. Clearly, this method probes 
the dark energy via both the growth rate and the z(r) relat ion. The method has recently been 



applied for the first time, to the COMBO- 17 survey (|l61|) . 
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4.6 Parameter forecasts from 3D lensing methods 



In this Section we summarise some of the forecasts for cosmological parameter estimation 
from 3D weak lensing. We concentrate on the statistical errors which should be achievable 
with the shear ratio test and with the 3D power spectrum techniques. Tomography should 
be similar to the latter. We show results from 3D weak lensing alone, as well as in com- 
bination with other experiments. These include CMB, supernova and baryon oscillation 
studies. The methods generally differ in the parameters which they constrain well, but also 
in terms of the degeneracies inherent in the techniques. Using more than one technique can 
be very effective at lifting the degeneracies, and very accurate determinations of cosmo- 
logical parameters, in particular dark energy properties, may be achievable with 3D cosmic 
shear surveys covering thousands of square degrees of sky to median source redshifts of 
order unity. 

The figures [8] and |9] show the accuracy which might be achieved with a number of surveys 
designed to measure cosmological parameters. We concentrate here on the capabilities of 
each method, and the methods in combination, to constrain the dark energy equation of 



state, and its evolution, parametrised by (|46|) 



w{a) = ^=WQ + Wa{\-a) (4.20) 

where the behaviour as a function of scale factor a is, in the absence of a compelling theory, 
assumed to have this simple form. The constant value w = — 1 would arise if the dark energy 
behaviour was actually a cosmological constant. 

The assumed experiments are: a 5-band 3D weak lensing survey, analysed either with the 
shear ratio test, or with the spectral method, covering 10,000 square degrees to a median 
redshift of 0.7, similar to the capabilities of a groundbased 4m-class survey with a several 
square degree field; the Planck CMB experiment (14-month mission); a spectroscopic sur- 
vey to measure baryon oscillations (BAO) in the galaxy matter power spectrum, assuming 
constant bias, and covering 2000 square degrees to a median depth of unity, and a smaller 
z = 3 survey of 300 square degrees, similar to WFMOS capabilities on Subaru; a survey of 
2000 Type la supernovae to z = 1.5, similar to SNAP'S design capabilities. 

We see that the experiments in combination are much more powerful than individually, as 
some of the degeneracies are lifted. Note that the combined experiments appear to have 
rather smaller error bars than is suggested by the single-experiment constraints. This is 
because the combined ellipse is the projection of the product of several multi-dimensional 
likelihood surfaces, which intersect in a small volume. (The projection of the intersection 
of two surfaces is not the same as the intersection of the projection of two surfaces). The 
figures show that errors of a few percent on wq are potentially achievable, or, with this 
parametrisation, an error of w at a 'pivot' redshift of z ~ 0.4 of under 0.02. This error is 
essentially the minor axis of the error ellipses. 
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Fig. 8. The accuracy expected from the combination of experiments dedicated to studying dark 
energy properties. The equation of state of dark energy is assumed to vary with scale factor a as 
w(a) = Wo + Wa{l — a), and the figures show the 1-sigma, 2-pai-ameter regions for the experiments 
individually and in combination. The supernova study fills the plot, the thin diagonal band is Planck, 
the near-vertical band is BAO, and the ellipse is the 3D lensing power spectrum method. The small 
ellipse is the expected accuracy from the combined experiments. From Heavens et al. (2006). 

4. 7 Intrinsic alignments 

The main signature of weak lensing is a small alignment of the images, at the level of a 
correlation of ellipticities of ~ 10^^. One might be concerned that physical processes might 
also induce an alignment of the galaxies themselves. The possible effect is immediately 
apparent if one considers that the shear is often estimated from the ellipticity of a galaxy, 
which includes the intrinsic ellipticity of the source Cs '■ 

e~e5 + 2Y. (4.21) 

A useful statistic to consider is the shear correlation function, which would normally be 
estimated from the ellipticity correlation function: 

(ee*) =4(Yy*) + (e.e:)+4(Ye*). (4.22) 
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Fig. 9. As in Fig. ID but with the shear ratio test as the lensing experiment. Supernovae fill the plot, 
Planck is the thin diagonal band, BAO the near-vertical band, and the shear ratio is the remaining 
45 degree band. The combination of all experiments is in the centre. From Taylor et al. (2006) 

This equation is schematic, referring either to galaxies separated by some angle on the sky, 
or by a 3D separation in the case of a 3D analysis. The first term is the cosmic signal one 
wishes to use; the second term is the intrinsic alignment signal, and the third is the shear- 
intrinsic alignment signal, which we will consider later. Until recently, both these additional 
terms were assumed to be zero. The hope was that even galaxies close together on the line 
of sight would typically be at such large physical separations that physical processes which 
could correlate the orientations would be absent. However, the lensing signal is very small, 
so the assumption that intrinsic alignment effects are sufficiently small needs to be tested. 
For the intrinsic alignment si gnal , this was first done in a series of papers by a number of 



t si gna 
59l;l42 



groups in 2000-1 (|96l : 16 11 : 1591 : 1421) . and the answer is that the effect may not be negligible, 
and is expected to be strongly dependent on the depth of the survey. This is easy to see, 
since at fixed angular separation, galaxies in a shallow survey will be physically closer 
together in space, and hence more likely to experience tidal interactions which might align 
the galaxies. In addition to this, the shallower the survey, the smaller the lensing signal. 
In a pioneering study, the alignments of nearby galaxies in the SuperCOSMOS survey 



were investigated (|38|) . This survey is so shallow (median redshift ~ 0. 1) that the expected 
lensing signal is tiny. A non-zero alignment was found, which agrees with at least so me o f 
the theoretical estimates of the effect. The main exception is the numerical study of (TO), 
which predicts a contamination so high that it could dominate even deep surveys. For deep 
surveys, the consensus is that the effect is expected to be rather small, but if one wants to 
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use weak lensing as a probe of subtle effects such as the effects of altering the equation of 
state of dark energy, then one cannot ignore it. There are essentially two options - either 
one tries to calculate the intrinsic alignment signal and subtract it, or one tries to remove 
it altogether. The former approach is not practical, as, although there is some agreement as 
to the general level of the contamination, the details are not accurately enough known. The 
latter approach is becoming possible, as lensing surveys are now obtaining estimates of the 
distance to each galaxy, via photometric redshifts (spectroscopic redshifts are difficult to 
obtain, because one needs a rather deep sample, with median redshift at least 0.6 or so, 
and large numbers, to reduce shot noise due to the random orientations of ellipticities). 
With photometric redshifts, one can downweight or completely rem ove physi cally close 
galaxies from the pair statistics (such as the shear correlation function) (|100l : ll59|) . Thus one 
remo ves a systematic error in favour of a slightly increased statistical error. The analysis in 
(|l01|) explicitly removed close pairs and shows that it can be done very successfully. 



4.8 Shear-Intrinsic alignment correlation 



The cross term {ye*) was neglected entirely until it was pointed out dlOSh that it was not 
necessarily zero. The idea here is that the local tidal gravitational field contributes to the 
shear of background images, and if it also influenced the orientation of a galaxy locally, then 
it could induce correlations between foreground galaxies and background galaxies, even 
though they may be physically separated by gigaparsecs. This term is more problematic for 
cosmic shear studies, because it is not amenable to the simple solutions which work well 
for the intrinsic alignment signal. It is conceivable that this effect is the limiting systematic 
effect in cosmic shear studies, as it seems necessary actually to model it and remove it. 
Studies of the SDSS (.179.) measured a significant signal in a related statistic, for very 
luminous galaxies, and a study of N-body simulations supported the view that the effect was 
likely to be non-negligi ble, at the l evel of up to 10% of the cosmic shear signal. On a more 
positive note, it seems dlOSl : II6OI : Il03|) that the term scales with source and lens angular 
diameter distances in proportion to the lensing efficiency (D (%/) © {Xs — Xi)/l> (%.?)■ This 
is reasonable, and also very useful, as it makes the parametrisation of the shear-intrinsic 
alignment much more straightforward. One can either use templates (|160|) or parametrise 
the contamination as a single function of separation, and marginalise over these nuisance 
parameters in the estimation of cosmological parameters. 
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5 Non-Gaussianities 



The two-point statistics discussed in Section [3] can be used to constrain cosmological pa- 
rameters. However, since they can be expressed in terms of the convergence power /k(^) 
they mainly depend on the same combination of parameters. Thus, from Eg. (12.91 ) we can 
expect (ic^) ~ o^Q.^ if we neglect the dependence on cosmology of comoving distances, 
where Os is the normalization of the linear power- spectrum. A more careful analysis (|22|) 
actually gives the scaling of Eq. (|10.1l) . In order to lift this degeneracy between the pa- 
rameters and Og one can combine weak lensing observations with other cosmological 
probes such as the CMB, as we shall discuss in Section [lOl or use 3D information as seen 
in Section|4](e.g. Eg. (14. 161 )). An alternative procedure is to consider higher-order moments 
of weak lensing observables. Indeed, even if the initial conditions are Gaussian, since the 
dynamics is non-linear non-Gaussianities develop and in the non-linear regime the density 
field becomes strongly non-Gaussian (this is an unstable self-gravitating expanding sys- 
tem). This can be seen from the constraints (5) = and 5 > — 1 (because the matter density 
p is positive) which imply that in the highly non-linear regime ((5^) ^ 1) the probability 
distribution of the density contrast 6 must be far from Gaussian. Since weak gravitational 
lensing effects arise from the matter distribution (see Eg. (12. 41) ) high-order correlation func- 
tions of both the 3D density field and weak-lensing observables are non-zero and could be 
used to extract additional information. 



5.1 Bispectrum and three-point functions 



The three-point correlation function is the lowest-order statistics which can be used to 
detect non-Gaussianity. In Fourier space it is called the bispectrum which is defined as: 

(5(ki)5(k2)5(k3)) = (27r)35z)(ki+k2+k3)5(fci,fc2,fc3) (5.1) 

for the 3D matter density contrast, where the Dirac factor results from statistical homo- 
geneity. Isotropy also implies that 5(^1,^27^3) only depends on the length of the three 
wavenumbers ki,k2,k3 or alternatively on two lengths ^1,^2 and the angle ai2 between 
both vectors. A key feature of the bispectrum (15.11) is that in the large-scale limit its de- 



pendence on the normalization of the power- spectrum can be factorized out (|22l) . Indeed, 
at large scales where the density contrast is much smaller than unity and guasi-linear per- 
turbation theory is valid one can expand the density contrast as a perturbative series of the 
form: 

5(k,z) = 5(i)(k,z) + 5(2)(k,z) + ... (5.2) 

where 5*^^) is of order q over the initial density field (6*^^-' is simply the linear density contrast 
6l). Then, substituting into the three-point function we obtain: 

(5(ki)5(k2)6(k3)) = (5(i)(ki)5(i)(k2)6(i)(k3)) + (5(2)(ki)6(i)(k2)6(i)(k3)) 
+(5(i)(ki)5(2)(k2)5(i)(k3) + (5(i)(ki)5(i)(k2)5(2)(k3) + ... (5.3) 
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where the dots stand for terms of order (5^^^)^ and beyond. For Gaussian initial conditions 
the first term vanishes whereas the three other terms are of order (S*^^^)^ so that the quantity: 

P{k2)P{h)+P{h)P{h)+P{h)P{k2) ^ ^ ^ 

is independent of the normalization of the linear density power- spectrum Piik) at large 
scales. In this manner one can separate the dependence on Cg from the dependence on 
other cosmological parameters. Using the small-angle approximation the £-space three- 
point correlation of the convergence reads (|27b : 

(k(/i)k(£2)k(£3)) = {2%)^Wi +h + i3)BAiuh,i3), (5.5) 

with: 



B^{h,h,h)= /dX-45 . (5.6) 



©4 X-D'CD'CD 



Then, as in Eg. (15.41 ) one can consider ratios such as 5k(|i,^2,^3)/(^k(^2)A:(^3) + ••) to lift 



the degeneracy between the parameters Qm and Og (|22D . Using tomography (i.e. redshift 
binning of the sources) also helps to constrain cosmological parameters such as the equation 
of state of the dark energy component, as studied in Ref. (i276) . We display their results 
in Fig. \W\ which shows that bispectrum tomography can improve parameter constraints 
significantly, typically by a factor of three, compared to just power spectrum tomography. 

In practice, most of the angular range probed by weak lensing surveys is actually in the tran- 
sition domain from the linear to highly non-linear regimes (from 10' down to l'). Therefore, 
it is important to have a reliable prediction for these mildly and highly non-linear scales, 
once the cosmology and the initial conditions are specified. Since there is no rigorous an- 
alytical framework to fully describe this regime numerical simulations play a key role to 
obta in th e non-linear evolution of the matter power spectrum and of higher-order statis- 
tics dlo^; I266). Based on these simulation results and analytical insight it is possible to 



build analytical models which c an de scribe the low order moments of weak lensing ob- 



servab les such as the bispectrum (12961) . Using a halo model as described in Appendix I A. 3 [ 
Refs. davi; |275|) investigated the real- space three-point correlation of the convergence 
(k(0i)k(02)k(03)). They studied its dependence on the triangle geometry (0i,02,93) and 
on the parameters of the halo model ( 274) and compared these predictions with numerical 
simulations (275,). 

As seen earlier it is more convenient for observational purposes to consider the shear rather 
than the convergence since it is the former which is directly measured (in fact what is 
actually measured is t he re duced shear y/(1 — k) which can be approximated by y in the 



weak-lensing regime, (|233[ )). However, since y is a 2-component field there are many ways 
to combine shear triplets. Here we defined the shear spin-2 "vector" as y = yiex + Ji^y 
where ex,ey are the 2D basis vectors using the fiat sky approximation which is valid for 
small angles (let us recall that y is not truly a vector since its components change as cos(2\|/) 
and sin(2\|/) under a rotation of \\r of coordinate axis, as seen from Eg. (12. 171) ). Besides, one 
must take care not to define statistics which depend on the choice of the coordinate system. 
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Fig. 10. Projected 68% confidence level constraints in the parameter space of Hje, wq, Wa and 
08 from the lensing power spectrum and the bispectrum in two redshift bins, as indicated. The 
coefficients wq and Wa parameterize the equation of state of the dark energy component. The results 
shown are obtained assuming priors on n, £lbh^ and h expected from the Planck mission. The sky 
coverage and number density are taken to be /sky = 0.1 and ?ig = 100 arcmin^^, and angular modes 
50 < / < 3000 are used. From Takada & Jain (2004). 

A possible approach is to consider scalar quantities such as the aperture mass M^p which 
can be expressed both in terms of the convergence or shear fields. However, this may not 
be optimal from a signal-to-noise perspective since the integration over the window radius 
Bs may dilute the cosmological signal as contributions from triangle configurations where 
the shear three-point function is positive or negative can partly cancel out. Moreover, the 
additional information contained in the detailed angular behavior of the shear three-point 
correlation can be useful to constrain cosmology and large-scale structures. Therefore, it is 
interesting to build estimators designed for the high-order correlations of the shear field. 

One strategy investigated in Ref.(|27h is to study the mean shear pattern y(0) around a pair 
of points 01,92 through the quantity {J3Q)) = ([y(0i)-Y(02)]y(0))- This study, based on an- 
alytical results (using the behavior of the density three-point correlation in the quasi-linear 
regime and its simplest extension to smaller scales) and numerical simulations, shows that 
this mean shear is almost uniform, and perpendicular to 612, over an elliptic area that covers 
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Fig. 11. Results for the VIRMOS-DESCART survey for the two point correlation function (left) and 
the reduced three point function (right). The solid line with error bars shows the raw results, when 
both the E and B contributions to the two-point correlation functions are included. The dot-dashed 
line with error bars corresponds to measurements where the contribution of the B mode has been 
subtracted out from the two-point correlation function (but not from ^3 there is no known way to do 
it). These measurements are compared to results obtained in xCDM, OCDM and ACDM simulations 
(dashed, dotted and dot-dashed lines respectively). From Bernardeau et al. (2002). 

the segment 0i2 which joins both points. This suggests to measure the average of (73(6)) 
over this ellipse so as to avoid cancellations. In this manner (2^ managed to obtain from 
the VIRMOS-DESCART Lensing Survey the first detection of non-Gaussianities in a weak 
lensing survey. We display their results in Fig. [TTI which shows that the amplitude and shape 
of the signal agree with theoretical predictions from numerical simulations. Although the 
measures are still too noisy to provide useful constraints on cosmology they show such 
weak-lensing observations to be a very promising tool. On the other hand, (28) also ob- 
tained explicit analytical expressions for the shear three-point correlations from the one- 
halo term which appears within halo models (when all points are assumed to lie within the 
same dark matter halo, see Appendix IA.3I) and recovered the pattern shown by numerical 
simulations. These results may serve as a guideline to build optimized estimators for the 
shear three-point correlations. 



A more systematic approach presented in (l243h is to look for natural components which 
transform in a simple way through rotations. Thus, to handle the three point function 
(Y(0i)Y(02)y(93)) one first defines the "center" c of the three directions 0i, 02, 63 on the sky 
(c may be taken for instance as the centroid, the circumcenter or the orthocenter of the tri- 
angle). Next, at each point one defines as the component of the shear along the direction 
that separates c and 0/ and Yx as the component along this direction rotated by 45°. From 
these tangential and cross components one introduces the complex shear y'^'^^ = Y+ + 'Yx 
and the "natural components" are defined as the four complex combinations: 



r(2) 



(y(c)(0j)y(c)(e2)Y(^)(03)), r(i) 



^y(c)*y(c)yc)^^ 



f^)yiC-)*y{c)^^ p(3) ^ 



c)y(c)y(c)*^_ 



(5.7) 
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Clearly the T^'^ only depend on the geometry of the triangle but to avoid ambiguities and 
miscalculations one must ensure that the points are always labeled in the same direction 
(e.g. counterclockwise, (244)). Each of these T^'^ is invariant only under special rotations 
but the important feature is that, under a general rotation, the different T^'^ do not mix but 
are simply multiplied by a phase factor. Note however that the four T^'^ are not independent 
as they arise from the same matter d istrib ution and three-point statistics are fully described 
by the projected matter bispectrum (|245h . On the other hand, such interrelations provide a 
redundancy which might be used to detect noise sources or B modes. 



An alternative method to study the shear three-point function is to divid e the possible com- 
binations into even and odd quantities through parity transformations (|312h . For instance, 
one chooses for the center of the triangle the barycenter c = (0i + 02 -l- 03)/3 and defines 
again the tangential and cross components of the shear, Y+ and Yx , from the direction that 
separates c and 0/ and from this direction rotated by 45°. Clearly this rotation changes 
direction under a parity transformation so that Yx — Yx (since Y is a spin-2 field, see 
Eg. (12. 171) . and the relative rotation is 2 x 45° = 90°) whereas Y+ Y+- Therefore, we ob- 
tain four parity-even three-point correlations: (Y+Y+Y+)' (Y+YxYx), (YxY+Yx), (YxYxY+)' 
and four parity-odd three-point correlations: (YxYxYx), (YxY+Y+)' (Y+YxY+)' (Y+Y+Y x ) ■ A s 
a consequence, for some symmetric configurations some odd functions must vanish (12731) . 
In particular, for equilateral triangles all odd functions vanish. This property assumes that 
the shear results only from weak-lensing (which only produces E modes), whereas source 
galaxy clustering, intrinsic alignments and observational noise can produce both E and B 
modes (see Sections [3] and |4l). The advantage of this procedure is that by focusing on even 
functions one avoids to dilute the signal by combining the estimators with parts which con- 
tain no weak-lensing information (odd functions for symmetrical triangle geometries). Be- 
sides, the parity-odd functions can be used to monitor the noise or to estimate the contribu- 
tion associated with higher-order effects beyond the Born approximation, source clustering 
or intrinsic alignments. Ref. (312) used a halo model (see Appendix |A. 31 ) to investiga te the 
behavior of these three-point correlations as a function of the triangle geometry. Ref. (l275h 
found that the halo model agrees well with numerical simulations at scales > 1' and could 
be used to obtain predictions for shear statistics in order to lift degeneracies in cosmolog- 
ical parameters. We display their results in Fig. [T2] which also shows that odd functions 
are smaller than even ones and vanish for symmetric geometries. They also note that fu- 
ture weak-lensing observations may be able to constrain the parameters of the halo model 
such as the mean halo den sity p rofile and halo mass function. On the other hand, using ray- 
tracing simulations Ref. (127 3h also evaluated the signal-to-noise taking into account the 
noise associated with galaxy intrinsic ellipticities. They found that a deep lensing survey of 
area 10 deg^ should be sufficient to detect a non-zero signal but an accurate measure would 
require an area exceeding 100 deg^. 
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Fig. 12. The eight shear 3-point correlation functions for the ACDM model against triangle config- 
urations (\|/ is the angle between the two sides of length q and r). The upper and lower plots show 
the results for the parity-even and -odd functions, respectively. Note that range on the y-axis for the 
right panel is about two times smaller than in the left panel. The solid curves show the halo model 
predictions for the eight shear 3-point correlation functions, while the symbols are the simulation 
results as indicated. From Takada & Jain (2003). 

5.2 Cumulants and probability distributions 



A simpler quantity than the three-point functions discussed above is provided by the third- 
order cumulant {X^) of smoothed weak-lensing observables (12.1 II) such as the smoothed 
convergence ic or the aperture-mass Map. In the quasi-linear regime where a perturbative 
approach is valid (with Gaussian initial conditions) one can see from Eqs. (|5.3l) - (I5.6I) that 

the skewness S^^^ = (k^)/ (ic^)^ of the smoothed convergence is independent of the mat- 
ter density power spectrum normalization as (the same property is clearly valid for other 
observables like Map which are linear over the matter density field). Therefore, by mea- 
suring the second- and third-order moments of the convergence or of the aperture mass 
at large angular scales one can obtain a constraint on O.^^ (|22h . Indeed, from Eq. (|2.9l) we 
see that k ~ Q.^ (neglecting the dependence of cosmological distances © on f^m) hence 
we can expect a strong dependence on Q.^ of the skewness as ~ • ^ numeri- 
cal study shows indeed that for sources at redshift Zs — 1 the skewness scales roughly 
as Sg^"* ~ ^^m*^'^ It^. Alternatively, from the skewness of weak lensing observables one 
can derive the skewness of the matter density field in the linear regime and check that 
the scenario of the growth of large-scale structures through gravitational instability from 
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initial Gaussian conditions is valid. In order to increase the information content which 
can be extracted from low-order cumulants one can consider generalized moments such as 
(A/ap(0.vi)A/ap(052)^^ap(0.s3)) which cross-correlate the aperture-mass Map(0.v,) associated 
with three different filter radii G^i, 0.v2 and 0^3- Then, the amplitude of such cumulants can 
be used to constrain cosmological parameters whereas the dependence on the angular ra- 
dius or the angular sepa ration of th e various filters helps constraining the properties of the 



large-scale density field (|203Hl58|) . For instance, one can introduce correlation coefficients 
rpqs such as: 

^^■^ ^j^p+q+sy^J {p+q+s) ^j^p+q+s\^q/{p+q+s) ^j^p+q+sy/ (p+q+s) 

where Xj is the aperture-mass or the convergence smoothed over scale 0„ (the source red- 
shift distributions ni{zs) may also be different). These correlation coefficients describe the 
information associated with three-point cumulants {X^X2X^)c which goes beyond the one- 
point cumulants {XP)c. We show in Fig. [13] the predictions of an analytical model based 
on a hierarchical ansatz (Appendix I A. 21) for the aperture-mass statistics, applied to the 
planned SNAP survey (left panel) and compared with numerical simulations (right panel). 
The behavior of these correlation coefficients can be used to discriminate between mod- 



els of the density field (12031) and to check that the observed non-Gaussianities arise from 
non-linear gravitational clustering. On the other hand, the left panel in Fig. [13] also shows 
how the aperture-mass is correlated between different angular scales. This correlation de- 
creases faster than for the smoothed convergence or the smoothed shear because the filter 
^Map(^0i ) is more narrow than W-^ in Fourier space, see Fig. [2] This is actually useful if one 
intends to derive constraints on cosmology from weak lensing surveys since it means that 
the errors associated with sufficiently different scales are uncorrelated. This also holds for 
the two-point moments discussed in Section[3]and for higher-order moments. Of course, the 
power- spectrum /k(^) and higher-order generalization such as the bispectrum5K(£i, ^2, ^3) 
are even less correlated and contain all the relevant information. 

In practice, because of the numerous holes within the survey area one first computes shear 
three-point correlations by summing over galaxy triplets and next writes (M^p) as an inte- 
gral over these three-point correlations (using the fact that Map can be written in terms of 
the shear y), see for instance (144). Applying this method to the VIRMOS-DESCART data 



Ref. (|213h were able to detect sf^"'^' and to infer an upper bound Q.^ < 0.5 by comparison 



with simulations. Next, one could measure higher order moments of weak lensing observ- 
ables. Note that for the shear components odd order moments vanish by symmetry so that 
one needs to consider the fourth-order momen t to g o beyond the variance (273). However, 



higher order moments are increasingly noisy (|291|) so that it has not been possible to go 
beyond the skewness yet. 

In order to compare observations with theory one needs to use numerical simulations or to 
build analytical models which can describe the low order moments of weak lensing observ- 
ables or their full probability distribution, as described in Appendix [A] This can be done 
through a hierarchical ansatz where all higher-order density correlations are expressed in 
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Fig. 13. Left panel: the correlator rm of the aperture-mass for the full wide SNAP survey is plotted 
as a function of smoothing angle 6.53 for a fixed pair of (6.«i,6i2)- The pair (641,6.42) for each curve 
is indicated in the plot. Error bars denote the 1 — a scatter around the mean, associated with galaxy 
intrinsic ellipticities and cosmic variance. Right panel: the three-point correlation coefficient ry\ \ of 
the aperture-mass as a function of smoothing angle 643 for a fixed pair of (6^1 , Qsi)- The three source 
redshifts are equal: Zs\ = Zsi = 2*3 = 1- The solid curve is the analytical model dA. 191 ) while solid 
points with error-bars are measurements from simulation data. From Munshi & Valageas (2005b). 



terms of the two-point correlation (13; 29 ih. Then, the probability distribution of weak lens- 
ing observables can be directly written in terms of the probability distribution of the matter 
density (Appendix IA.2I) . In some cases the mere existence of this relationship allows one t o 
discriminate between analytical models for the density field which are very similar (|20l[) . 
Alternatively, one can use a halo model (Appendix IA.3I) where the matter distribution is 
described as a collection of halos and the low order moments of weak len sing ob- 
servables can be derived by averaging over the statistics of these halos ( 2741 : |275 ). On the 
other hand, one can use weak lensing to constrain halo properties and to detect substruc- 
tures (j66j). For the particular case of the smoothed convergence the probability distribution 
function (PDF) T (ic) can be expressed in terms of the PDF T (5) of the 3D matter density 
contrast within some simple approximations, see Eq. dA.VI) . This allows one to apply to the 
convergence simple mod els w hich were originally devised for th e 3D density field such 



as the lognormal model (|278h or more elaborate ones (|287l : l305l : Il3h . For more complex 



weak-lensing observables such as the aperture-mass these approximations can no longer be 
used and one needs an explicit model of the density correlations (see AppendixjAj) to derive 
their cumulants or their PDF We display in Fig.[l4]the results obtained from a hierarchical 
ansatz f or th e aperture-mass, for the skewness (left panel) and the PDF of an estimator 
of Map (|29lh . In such calculations one must take into account the noise associated with 
the intrinsic ellipticities of galaxies and the cosmic variance. However, the left panel of 
Fig. [m shows that despite these sources of noise future surveys such as SNAP should be 
able to constrain cosmological parameters at a level of 10% from such low-order moments, 
whereas the right panel shows that the tails of the PDF '2 (Map) should allow one to extract 
information beyond low-order moments. 
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Fig. 14. Left panel: The skewness 5^"'' = (M^p) / (M^p)^ of the aperture-mass (sohd curve), for the 
SNAP survey. The central error bars show the 1 — a dispersion due to galaxy intrinsic ellipticities 
and cosmic variance. The smaller error bars which are slightly shifted to the left show the dispersion 
obtained by neglecting non-Gaussian contributions to the dispersion whereas the smaller error bars 
which are slightly shifted to the right show the dispersion obtained from the estimator built from the 
cumulant rather than the moment. We also show the effect of a 10% increase of Q.^ (lower dotted 
curve), of a 10% increase of ag (central dot-dashed curve) and of a 10% decrease of the redshift zq 
(upper dashed curve). Right panel: The pdf 9 (M) for the estimator M associated with the aperture- 
mass Map. Note that the Gaussian noise introduced by intrinsic ellipticities makes !p(M) closer to 
the Gaussian than the actual pdf !p(Map) which only takes into account gravitational lensing. The 
solid line shows the theoretical prediction, the dashed line is the Gaussian and the dotted line is 
the Edgeworth expansion up to the first non-Gaussian term (the skewness). The error bars show the 
1 — a dispersion. From Valageas, Munshi & Barber (2005). 

5.3 Primordial non-Gaussianities 



So far we have discussed the non-Gaussianities associated with the non-linearity of the 
gravitational dynamics, assuming Gaussian initial conditions. However, results from future 
surveys can also be very use ful in constraining primordial non-Gaussianity predicted by 
some early universe theories On the other hand, generalised theories of gravity can 

have very different predictions regarding gravity-induced non-Gaussianities as compared 
to General Relativity, which can also be probed using future data. A joint analysis of power 
spectrum and bispectrum from weak lensing surveys will provide a very powerful way to 
constrain, not only cosmological parameters, but early universe theories and alternative 
theories of gravitation, see also (229). Thus, a recent work (28) has studied the possibil- 
ity of constraining higher-dimensional gravity from cosmic shear three-point correlation 
function. 
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6 Data Reduction from Weak Lensing Surveys 



6.1 Shape measurement 



Weak lensing by large scale structures induces a coherent alignment of galaxy shapes across 
large angular distances. A crucial step for its measurement is the accurate estimation of 
galaxy shapes, free of biases and systematics. A common approximation is to describe 
galaxies as simple elliptical objects, for which the quadrupole of the light distribution is a 
fair estimate of the shear y. If we call /(0) the 2-dimensional light distribution of the galaxy 
image, then the quadrupole moment Qij is defined in its simplest form as: 

_ /d9,d9,/(9)(9,-9,)(9,-9,) 
/ d9,d9, /(9) 

where 9/ is the centroid of the light distribution. The ellipticity e,- = (ei,e2) of the galaxy 
is given by (using the same definition as in Section [34l with Eqs. (l3.15l) - (|3.17l) ): 

^ ' Veil +322 211+622/ 

When galaxies are not lensed (in the source plane), their average ellipticity is zero (e,) = 
from the assumption of isotropy. Lensed galaxies exhibit an average non-zero distortion 
6/ = {et) over the coherence scale of the lensing potential (|) because the shear y,- stretches 
galaxy shapes locally in the same direction. The distortion is given by Eq. (l6.2l) . where 
Q is the quadrupole moment of the galaxy image . The relation between the source 
quadrupole and depends on the magnification matrix ^ (see Section O: 

<2^= <^-i<2^^-i, (6.3) 

where 

( 1-K-yi -y2 \ 
Si=\ . (6.4) 

y -y2 1 -K + yi J 

The matrix ^ of Eq. (l6.4l) describes the mapping between the source and image planes. 
The source galaxy surface brightness /■^(9) is stretched along the eigen axis of J^, and the 
observed galaxy light distribution becomes: 

/(9)= [l + (^-/),^.9;a,]/(9) (6.5) 

Figure [15] shows what is happening to the shape of a galaxy in the weak lensing regime. 
The quantities of interest can be listed as follows (see also Section O: 

• Convergence K 

• Distortion |6|=4tii 
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Fig. 15. Schematic representation of a lensed galaxy. In the source plane, the galaxy is circular and 
has radius Rq. Convergence K stretches its average radius to /?o/(l — k), and the shear y distort the 
galaxy along some angle (3. a and b are the semi-major and semi-minor axis respectively. 

. Shear |y|=^ 

• Magnification /i = [det(j?)]"^ = [(1 - k)^ - y^] 

The convergence and magnification cannot be easily measured because we do not know 
the size of galaxies in the source plane. Therefore, we focus on the shear (and the related 
aperture-mass Map, see Eg. (12. 121 )) as estimators of weak lensing effects. However we would 
like to point out that magnification has been measured successfully on SLOAN from the 
excess of distant quasars around foreground galaxies ( 250) and this technique could poten- 
tially be a powerful probe of cosmology in the future. The distortion offers the advantage 
that we know statistically its value in the source plane, i.e. (5,) = 0, meaning that source 
galaxies are randomly aligned (but see Sections 1477] and WM . 



6.2 Point Spread Function correction 



One cannot easily measure the distortion 5/ for two reasons: the optics of the telescope 
induce geometrical distortions that can be misinterpreted as a weak lensing signal if not 
carefully corrected: the stars, which provide a picture of the Point Spread Function (PSF) 
of the telescope's aperture, have complicated, elongated shapes. The second problem is 
that the weak lensing signal is better measured on distant galaxies. Indeed, since the line 
of sight is more extended weak lensing effects have a greater magnitude, see Eq. (|2.4l) . 
However, these distant galaxies are small and faint, and the sky noise for these objects is 
large. Consequently, a naive estimation of the quadrupole Qij from Eq. (|6.1| ) is essentially a 
measure of sky noise. 



Refs. (|150l : |173|) showed that one can overcome both problems by 1) reducing the noise at 
the edge of the galaxy using a gaussian filter, 2) expressing the effect of the PSF convolution 
analytically as a perturbation expansion, and use the first-order term to correct the galaxy 
shapes. The KSB method (150) can be summarized as follows: the observed ellipticity e"^''^ 
is the sum of three terms, the first is the intrinsic ellipticity ef' of the galaxy before lensing 
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and before convolution with the PSF. ef' is unobservable, but the isotropy of space implies 



0. The second term describes the galaxy shape response to an anisotropic PSF 
(which depends on a measurable quantity called the smear polarizability tensor P-"^), and 
the last term describes the response to the isotropic PSF (which depends on the "preseeing" 
shear polarizability P^, also measurable). Therefore the observed ellipticity becomes 



(6.6) 



where y, is the shear we want to measure and p* is the stellar ellipticity estimated at the 
galaxy position. P™\ p* and P^ can all be estimated from the image (see (fl5d) for the 
detai l s). S ome refinements should be included when shapes are measured on space data 
(11091 : |218|) . but they do not change the philosophy of the method, which remains a first- 
order perturbative approach. The estimation of p* at the galaxy position is done by first 
measuring p* on the stars and then interpolating its value assuming a second-order polyno- 
mial variation across the CCDs. This is a crucial step since an inaccurate model could lead 
to significant B modes in the signal (|298|) . fort unatel y var ious models have been proposed 
to account for non-polynomial variations (|114l) . Ref. has shown that we can perform 
a singular value decomposition method of the PSF variation between individual exposures 
in order to improve the correction. This would be a particularly useful approach for lensing 
surveys planning to observe the same part of the sky hundreds of times (LSST, ALPACA). 



The ultimate accuracy of galaxy shape correction is still a wide open question. Whether 
or not there is a fundamental limit in the measurement of galaxy shapes is a particularly 
critical issue for the design of future gravitational lensing surveys (see Section [63l) . KSB 
has been historically the first shape measurement method working and for that reason it 
has been tested intensively, but its accuracy is not expected to be better than 5-10% (|72l) . 
Clearly this is not enough for precision cosmology which seeks sub-percent precision on 
shape measurement. 

Many post-KSB techniques have been developed over the past five years in order to im- 
prove upon the original KSB approach, most of them are being intensively tested only now 
(|l02l : ll88|) . A quick summary of these methods follows. The r eader can obtain the details 
by looking at the original papers or the description given in /lO^). Ref. (11631) proposed 
to model the PSF and the galaxies by a sum of Gaussians of different widths. The pre- 
seeing galaxy shap e is recovered by a %^ minimisation between the galaxy model and the 
measured profile. (Il53|) extended the original KSB by properly modelling the PSF with 
a realistic kernel, dropping the assumption of a Gaussian profile. The PSF is then circu- 
larized prior to measuring the galaxy shapes (the circularization technique is also used by 
f265.)). A radically different approach c onsi s ts in projecting the shapes (galaxies and stars) 
on a basis of orthogonal functions (|29|; l220l : l22l[) . The effect of convolution and shear can 
be expressed analytically on the basis functions and the solution of the preseeing galaxy 
shape can be found by a straightforward matrix inversion. This technique turns out to be 
important not only for shape measurement, but also for simulating realistic galaxy profiles. 
It also offers in principle a total control of the different processes changing the shape of a 
galaxy (shear, amplification, PSF, etc.). Mathematically, a galaxy with profile /(0) can be 



52 



decomposed over a set of basis functions B(^„^ (0; [3) as 

m= E /(„,,„,)5(„,,„,)(e;|3), (6.7) 

(ni,n2) 

where p is a scaling parameter adjusted to the size of the galaxy we want to analyse, 
f [,11,112) = /n ^re called the shapelets coefficients, and could be any family of poly- 

nomial functions fulfilling our favorite recurrence relation (orthogonality, orthonormality, 
etc.). Although the description is given here for a cartesian coordinate system, the same 
formalisr n can be developed for any coordinate system better suited for galaxy shape anal- 



ysis (see (|186|) for a derivation using polar coordinates). The convolution of a lensed galaxy 
/(0) with a stellar profile ^(0) produces a galaxy with profile h{Q) such that its shapelets 
coefficients are: 

hn = ^ Cnml/m^l = ^ Pnmfm , (6.8) 
m,l m 

where C and P = Cg are matrices that can be measured on the data. The preseeing galaxy 
profile can be formally obtained from the matrix inversion f = P^^h. (|29|) developed a 
similar method but they also perform a circulariz ation of the PSF like in pnor to 



the measurement of the pre- seeing galaxy shap e. (11641) has implemented a version of the 



shapelet technique based on ideas developed in (|163l) where galaxies and PSF are decom- 
posed on a fixed set of simple profiles. Currently, there are six different shape measurement 
techniques, and several implementations of KSB, corresponding to more than a dozen of 
different pipelines. This clearly demonstrates the richness of the topic, but one should en- 
sure that they all lead to the sam e she ar measurement. A comparison of the performances of 



all these techniques is shown in (11061) . but the authors focussed on analytical galaxy profiles 
instead of real, noisy, profiles. They found that KSB was performing the best, which is sur- 
prising gi ven the number of approximations involved. The Shear TEsting Program STEP 



(102; 188:f ^has been setup in order to systematically perform intensive tests of shape mea- 



surement methods. Its goal is to test the different pipelines on simulated and real data sets 
and to find the ultimate limit of shape measurement from space and ground based images. 



(|102l) measured the difference between the measured shear y^^^^ and the true shear yj™® 
on a large number of simulated images with very different PSF simulating various optical 
defects of the telescope. The chosen parametrization was: 



yrneas _ yrue _ ^^yrue^ 2 ^ ^yme ^ ^ 

where q, m and c\ are measured from different realizations of image quality (PSF anisotropy 
and seeing size). Figure [16] shows the r.m.s. of these parameters. It shows the main result 
from STEP, indicating that shape measurement accuracy is within the few percent range, 
at least one order of magnitude above the required accuracy for future experiments such 
as SNAP and LSST. This figure also shows that KSB is performing remarkably well, in 
agreement with (|1060 . and the newer methods are potentially even more powerful. 



http://www.physics.ubc.ca/~heymans 
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Fig. 16. From Heymans et al. (2006), plot showing the calibration bias m, PSF residual and non- 
linearity (q) for different pipelines (see Heymans et al. 2006 for a complete Ust of references and 
detailed description of the methods). 

Among all the methods, the very attractive feature of the shapelets lies in the fact that each 
transformation experienced by a galaxy (PSF, shear, convergence, etc..) can be expressed 
as a simple set of operators acting linearly on the set of basis functions describing the 
galaxy. A general transformation of the galaxy shape is therefore a linear process which can 
formally be solved in one pass to provide the pre-seeing shape. The shapelets, or similar 
approaches (l29h . provide in principle a perfect description of the galaxy shapes. Is it the 
ultimate method with the best possible accuracy? In the weak lensing regime for instance, 
the quadrupole of the light distribution fully describes the shear y, a two-components spin- 
weight 2 object, and it is unnecessary to measure the details of the galaxy shape. In that 
case, the shapelets might appear like overkill, and limiting the number of basis functions is 
appropriate (164). The solution of the shape measurement problem is a trade-off between 
an accurate description of the galaxy morphology (i.e. galaxy structures) and an unbiased 
measure of the second order moments: we want to describe the galaxy shape with enough, 
but not too many, details. This optimal trade-off depends on the weak lensing information 
we want to extract from the galaxy distortion. For instance, in the weak lensing regime, 
it is assumed that the shear does not vary across the galaxy, and the quadrupole of the 
light distribution is then a complete description of the lensing effect. This is equivalent 
to saying that the centroid of a lensed galaxy is the lensed centroid of the source galaxy. 
Mathematically, this means that 
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(6.10) 



where 0"^ and 0^ are the angular position of a galaxy source and image respectively, and 
0"^^ and 0^^ are the centroid position of the source and image respectively. Eq. (l6.3l) is a 
direct consequence of this approximation. It obviously breaks down near the critical line 
where the determinant of the amplification matrix Si is zero and the magnification is infi- 
nite. The source galaxy is then strongly distorted and gravitational arcs are observed (,75.) . 
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In the intermediate regime where source galaxies are mildly distorted and have an arc -like 
shape, weak lensing is just an approximation, and the quadrupole of the light distribution 
is not enough to quantify the lensing e ffect . A higher-order description of the galaxy shape 



becomes necessary (186(). Refs.dSTI: l9l: 11891) developed the theory of flexion which is a de- 



scription of the next order of shear measurements, the octopole. In that case, Eq. (|6.10l) is 
not valid and should be replaced by: 



e; = ^„0, + -A7^0,0;, (6.11) 



where Djjk = ^k^ij is given by: 
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Higher orders of g alaxy shapes are a probe of higher order derivatives of the gravitational 
lensing potential (|239|) . This is therefore particularly relevant for lensing by cluster of 
galaxies and space quality images, because the latter is well suited for an accurate mea- 
surement of galaxy shapes beyond the quadrupole. The new relation between the source 
and image galaxy profile is given by 
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(6.13) 



which replaces Eq. (l6.5l) . The flexion terms can be conveniently expressed in terms of 
shapelet operators (|87i) . The practical utility of flexion has not been demonstrated yet, but 
progress is being made to show whether or not it is measurable. With the development of 
CCD detectors in space it is likely that flexion could provide useful lensing information 



Lots of progress has been made in the measurement of galaxy shape s since the original KSB 
paper in 1995. The most recent shape measurement method shows that the one per- 
cent accuracy goal has not yet been reached, and most of the effort now consists of showing 
that this goal can be met in order to perform high precision cosmology with weak lensing; 
this is the primary goal of the STEP project. With the development of these techniques and 
improving image quality, it becomes possible to measure higher order shear effects such 
as the octopole. Shapelets or a similar method are particularly useful to extract high-order 
morpholog y information, while KSB is enough for percent precision shear measurement 
(fl02l : [im . 



The shapelets can also be used to generate realistic images of galaxies, which is important 
for testing shape measurement methods. Here one uses a training set to teach a galaxy 
image simulator how to generate a realistic combination of shapelet coefficients in order 
to reproduce statistically real data sets. One such galaxy simulation is shown in Figure [TTI 
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Fig. 17. Simulating galaxies with shapelets: top-left is a galaxy from the Hubble Ultra Deep Field 
and the other panels show its simulated image using the shapelets with nmax = 5, 10, 30. 

for instance simulated galaxies using the shapelets, whose distribution has been 
trained using the Hubble Ultra Deep Field. 



6.3 Statistical and Systematic Errors 



In this Section, we consider the potential sources of error and systematics in the shear 
measurement, not the error caused by intrinsic alignment (see Section 14.71) and selection 
biases correlated with the shear orientation and amplitude dlOq) . The latter was shown to 
be negligible (.102) . 



A complete description of galax y shape is certainly not necessary for most weak lensing 
applications. According to (11021) . the main limitation in shear measurement comes from 
the calibration of the shear amplitude: the PSF anisotropy is relatively easy to correct, but 
the isotropic correction due to the seeing is still not accurate to better than a few percent 
(see Figu re [T6l). An additive error should also be considered, as suggested by Figure [T6l 
Ref . (11321) has shown that, in order not to degrade significantly the cosmological parameters 
estimation from future weak lensing experiments (SNAP and LSST), the additive error has 
to be less then 10^^, which is one order of magnitude better than what can be achieved 
today (see Figure [T6l). It is still unclear how to estimate the impact of the additive error for 
realistic surveys, since the redshift, color and morphology dependence might be quite com- 
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Fig. 18. Two sources of error: statistical noise and cosmic variance. The curves show the diagonal 
part of the covariance matrix for either source of noise, normalised by the fiducial model shear 
variance. The thick solid lines show various levels of multiplicative errors from 0.01 to 0.001. 

plicated and are not yet well understood. The multiplicative (i.e. calibration) error of the 
order of on e per cent does not seem to degrade dramatically the cosmological parameters 
constramts (432*), but this is because the cosmological constraints come from the largest 
scales, w here t he error budget is dominated by cosmic variance and not by the multiplica- 



tive error (|300f) . Figure [T8] shows that a large multiplicative error degrades significantly the 
weak lensing signal at angular scales less than 10 arcminutes, and has no effect at scales 
20 arcminutes and above. Therefore a complete scientific use of a lensing survey is also 
dependent on ou r ability t o reduce the multiplicative error, not only the additive error. With 
the STEP effort 1188.) and the large amount of ground and space based data available 
from current weak lensing surveys, the shape measurement issue will probably be solved 
within the next 2 or 3 years. The major source of error for weak lensing surve y wi ll then 



become our ability to estimate the photometric redshifts of the source galaxies (13001) . This 
is a powerful technique, but the redshift error is large and often there are multiple redshift 
sol ution s due to spectral features not covered by the set of filters. This is clearly shown 



in (|133r) . where photometr ic re dshift degeneracies are left with the 5 filters ugriz of the 



MEGACAM camera. Ref. (ll32|) has shown that the requirement of photometric redshift 
precision is tight if we want to achieve tomography: for instance the average redshift needs 
to be accurate to better than 1%. A large number of filters in the optical and near infrared 
coverage would then be necessary. Ref. (.161.) finds that the requirements are less severe for 
3D weak lensing (see Section |4]). 
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7 Simulations 



In order to obtain good signal-to-noise for estimates of cosmological parameters, it is nec- 
essary to probe many different scales, including small scales where linear or second-order 
perturbation theory is not valid. As we have no exact analytical description of matter clus- 
tering at small or intermediate scales, numerical N-body simulation techniques are em- 
ployed to study gravitational clustering in an expanding background. Numerical techniques 
typically use ray-tracing techniques through N-body simulations to study weak lensing of 
background sources. Other methods include line-of- sight integration of shear. Although 
often only limited by computational power, numerical techniques too depend on various 
approximations which can only be verified by consistency checks against analytical results. 



Simulating (strong) lensing by individual objects can provide valuabl e informati on regard- 
ing background cosmology through the statistics of arcs, see e.g. (12841 : |204|) for early 
studies. Simulating weak lensing surveys on the other hand probes inhomogeneous mat- 
ter distribution by large scale structure in the Universe. Early attempts to simulate weak 
gravitational lensing by inhomogeneous dark mat ter distribution was initiated by various 
authors in the early 1990s. These studies include (l23ll : CT : [TesI : Isl: [itI: [32|)^owever the 
most detailed studies in th is directio n started emerging by the late 1990s. (|304|) following 
up previous work done by (|302l : l303h presented a detailed analysis of lensing by large scale 
structure. On a slightly larger angular scale this study was extended and complemented later 
by the work done by (137) who constructed shear and convergence maps using ray-tracing 
simulations. With ever increasing computational power the typical size of the sky which 
can be simulated using ray tracing experiments has increased over the years and recent 
studies can now focus even on scales comparable to tens of degrees while still resolving 
smaller angular scales well. These numerical developments have also stimulated improved 
analytical modelling of weak lensing using perturbative techniques at larger angular scales 
and halo-based models or the hierarchical ansatz at small angular scales. 



7.1 Ray tracing 



In case of experiments involving ray-tracing simulations, one combines several large-scale 
boxes obtained from cosmological N-body simulations to build a large simulated volume 
from the observer up to the source plane at redshift Zs- Then, the dark-matter particle dis- 
tributions whithin each box are projected on successive two-dimensional planes up to the 
redshift of the source. Typically up to 30 such planes are employed to sample the matter 
density to a source redshift of ~ 1 and 10^ rays are propagated through the N-body data 
volume. The computation of the derivatives of the gravitational potential in each of these 
lens planes is performed using FFTs. This provides the shear tensor at each plane which 
allows to follow the deflection of the light rays from one plane to the next. The Jacobian 
matrix Si of the mapping from source to the image plane determines lensing observables 
such as the convergence k and the shear components Yi, Y2, see Eq. (|2.6l ). Such an algorithm 
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Fig. 19. The blue elongated disks represent observed images of background galaxies. The dark mat- 
ter filaments from numerical simulations which were used to simulate the survey are also plotted to 
show alignments of ellipticites of observed galaxies with underlying filamentary structures (Figure 
courtesy: Stephane Colombi and Yannick Mellier). 



is also known as the multiple-lens algorithm (see e.g. (|232|) for more detailed discussions). 



Various numerical artefacts that determine the resolution of a ray-tracing simulation include 
the spatial and mass resolution of the underlying N-body simulation through which the ray 
tracing experiments are being performed as well as the size of the grid which is used to 
compute the intermediate projected densities and the gravitational potential. The finite size 
of the simulation boxes on the other hand determines the largest angular scales to which 
we can reliably use the results from ray-tracing simulations. Depending on the size of the 
simulation box and source redshift one can typically construct a few degree square patches 
of the sky. To improve the statistics, the N-body simulation box is rotated and ray tracing 
experiments repeated, to generate additional weak lensing sky patches. 



Studies using ray tracing simulations were initiated by (Il37|) . They used a 256^ adaptive 
P^M simulation outputs from Virgo to perform the ray-tracing simulations. The lens plane 
grid used to compute the potential and its derivatives from the projected matter distribution 
had a resolution of 2048^. Depending on a specific cosmology these studies generated weak 
lensing maps of a few degrees, with resolution down to sub-arc minutes. Recent ray tracing 
simulation using multiplane techniques include the ones presented in (92) 



7.2 Line-of-sight integration 

Ref.dssj) pointed out that several problems can arise in weak lensing studies from the use of 
the multiple-plane approach described in the previous Section, especially when the sources 
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are distributed at relatively high redshifts. In particular, the projection of the matter distribu- 
tion onto successive 2D planes orthogonal to the mean line of sight clearly approximates all 
angular diameter distances by a constant within a given redshift interval. This may lead to 
significant errors if these intervals are too large. Thus, Ref.ltssI) introduced a 3D algorithm 
which computes the second derivatives of the gravitational potential 4>(x) within the full 3D 
volume using FFTs. Then, light rays are followed within the 3D volume and deflections are 
taken into account along the 3D grid using the local derivatives of the gravitational poten- 
tial and the local angular distances, which allows to derive the Jacobian matrix between 
source and observer planes (for more detailed discussions see (|58l : llOl : lUj)). The variance 
of the cosmic shear obtained from this algorithm was compared with an alytical pr edictions 
and other simulations in (12). More detailed comparisons appear in (|13u29ll : l201|) where a 
good agreement was found with a whole range of analytical predictions. 



It should be noted that ray-tracing simulations and line-of-sight integrations still remain 
costly options to simulate a reasonable portion of the sky with fairly low resolution. This 
limits the number of independent realizations which can be simulated. However, to probe 
the small angular scales where non-linear gravity has not been understood analytically yet, 
simulations are the only reliable option against which all analytical predictions are tested 
regularly. 
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8 Weak Lensing at other wavelengths 



Current weak lensing surveys mostly rely on statistical studies of ellipticities of background 
galaxies at optical wavelengths. However various authors have considered the possibility 
of weak lensing studies in other wavebands, both for individual radio or IR sources at high 
reds hift a nd for the fluctuations in the integrated diffuse emission from unresolved sources 
(|55l : l214l) . It was pointed out that future facilities at radio wavelengths will even start com- 
peting with space-based optical observations which are limited by their ability to resolve 
the shape of distant sources and their small field of view. Radio surveys by the proposed 
Square Kilometre Array (SKA) will be abl e to make huge progress in this direction by re- 
solving orders of magnitude more sources (|238|) . In addition to resolved individual sources 
it was also suggested that integrated diffuse emission from the first stars and protogalaxies 
at high redshift {zs — 15 — 30) as well as 21cm emission by neutral intergalactic medium 
can provide useful arenas for weak lensing studies. It was realized that such programmes 
could be very useful in bridging the gap between weak lensing surveys based on optical 
studies of nearby galaxies at a redshift of a few and weak lensing studies of the CMB at 
a redshift of Zs = 1 100. Nevertheless separating galactic contamination from cosmological 
signal remains a difficult task. 



8.1 Weak lensing studies in Radio and near IR 



Future radio facilities, such as the proposed skaE], will be superior to current radio ob- 
servatories by orders of magnitude, particularly in its field of view, and in sensitivity. In 
addition, the higher resolution achieved by SKA will place it in a much better position than 
the present generation of radio telescopes. Surveys using SKA will push radio astronomy 
in a position where the number density of radio sources will be comparable to that of the 
optical sky (where the number density of useable galaxies is ~ 30 arcmin up to a red- 
shift of order unity, for ground-based surveys). This is possible thanks to the fact that SKA 
will be able to observe 100 times fainter objects than currently achievable; the radio sky at 
present is almost literally empty. For effective weak lensing studies the number density n of 
sources as well as the mean redshift of sources (zs) should be as high as possible. However 
predictions about these faint radio sources and their {zs) and n is currently less certain. If 
the dominant population consists of normal or star-forming galaxies the average redshift 
would be roughly unity. Additional populations of sources can make the redshift distribu- 
tion less certain. Another complexity in the whole scenario is that we do not know to any 
accuracy the shape distribution (i.e. the value of Og) of the faint radio sources. Ideally one 
would hope for near-spherical sources with no intrinsic ellipticities acting as source objects, 
i.e. with Oe as low as possible. If the radio sources in SKA are dominated by core-jet type 
objects, then Og can be very high, whereas normal galaxies have relatively lower values. 
The field of view (FOV) for SKA will be large and clearly the PSF will be controllable; 



^ http://www.skatelescope.org/ 
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this makes it comparable or probably better than the present generation of optical surveys. 
As discussed above, seeing provides a fundamental limitation for optical telescopes: as 
they become fainter, the source galaxies tend to become smaller too. This limits the num- 
ber density of sources for which the ellipticity can be measured reliably. A situation where 
point sources dominate resolved source number counts is already present in some of the 
deepest space-based images available to date. These include the Hubble Deep Fields North 
and South and recent images from the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). Clearly the space-based images are much better compared to their 
ground-based counterparts, but the FOV of such observations is limited. Radio observations 
using future facilities could be highly productive as they could combine high measurable 
source density and large FOV. However, these conclusions will depend to some extent on 
the intrinsic properties of the sources, and further assume that sources are not too elongated 
and that the average redshift distribution of these sources is not too shallow. 



8.2 Possibility of 21cm weak tensing studies 



Usually background objects can be broken down to indivi dual s ources at opti cal and infra- 
red wavelengths, but the emission of 21cm radiation ( 2481 : 1 1771 : 12831: 1 1 341: l77h from neutral 
gas prior to re i oniz ation, and the CMB, provide examples of truly diffuse backgrounds 



( 211 



2711 : 12631 : 1126 ) 



At high redshifts, typically beyond current large-scale surveys, the unresolved point sources 
become more dominant. At a source redshift of z^, = 15 — 30, first generations of stars 
and protogalaxies start to appear according to current theory of galaxy formation. These 
point-like objects provide a perfect background for weak lensing studies due to the large 
distances that light rays need to travel from these objects to reach us. Though detection 
of these objects is beyond present observational technology, the spatial fluctuations in the 
integrated diffuse background emission can be potentially interesting for lensing studies if 
we assume the most optimistic emission models. 

The analytical formalism for weak lensing studies of diffuse backgrounds has recently 
been established (55). Borrowing techniques from studies of weak lensing of the CMB, 
this shows that a perturbative approach - typically employed to study CMB lensing - will 
still be valid for diffuse background studies. This is despite the fact that unlike the CMB 
there is considerable power at small angular scales due to the lack of damping tail. This 
indicates that a perturbative series has a comparatively slow convergence rate and a larger 
number of terms need to be included for realistic calculations. 



In contrast to CMB lensing studies (I310l : ll9l : l9d : ll25l : ll07|) . the weak lensing studies based 
on diffuse components suffer from the fact that the lensing modification to the power spec- 
trum is minor at arc-minute angular scales and the lensing information that one can extract 
from low-redshift diffuse background is significantly limited. This is related to the pres- 
ence of significant structure in CMB power spectra, as signified by the acoustics peaks. 



62 




10* 10« 



Fig. 20. In the left panel convergence power spectra for redshift z.s = 9 are plotted. The dotted line on 
the left corresponds to noise level expected from PAST/LOFAR where as the dotted line on the right 
corresponds to noise level expected from second generation PAST+ and the SKA. Measurement 
errors and accuracy of various surveys in estimating matter power spectra are compared on the right 
panel. Solid line corresponds to the matter power spectrum at z=9 as in left figure. Dot-dashed line 
corresponds to power spectrum of patchy reionization with a bias factor b = 4 for more details). 
Dashed lines correspond to noise levels expected for PAST and second generation PAST+. Dotted 
line corresponds to fractional accuracy of power spectrum estimation per logarithmic I bins. (Figure 
Courtesy Ue Li-Pen) 

which is lacking in the case of diffuse background studies. The presence of a damping 
tail in the CMB means that the convolution associated with weak lensing effects transfers 
some power from larger angular scales to smaller angular scales at the arc-minute range, 
which makes it most easily detectable. Use of polarization information in CMB also carries 
much richer lensing information as compared to unpolarized diffuse background studies. 
Moreover, at the last scattering surface where the CMB temperature and polarization fluc- 
tuations are generated they follow a Gaussian pattern. Then lensing due to the intervening 
mass distribution imprints a no n-Ga ussian footprint which can be effectively used to extract 
cosmological information (23; 15 vl). There have been extensive studies in this direction us- 
ing higher-order moment-based techniques (see e.g. (54)). By contrast the non-Gaussianity 
generated by weak lensing in diffuse backgrounds such as the 21cm background is unlikely 
to be detected in near future. Nevertheless if possible 21cm weak lensing studies can extend 
the reconstruction of the integrated matter power spectrum out to redshifts of 15 to 30, and 
will bridge the gap between current and upcoming galaxy lensing studies and the CMB. 



8.3 Using resolved mini-halos for weak lensing studies 



Concentrating on mini-halos that will contribute to the 21cm background radiation, (^2141) 
found that instruments such as PAST, LOEAR would tremendously improve our knowledge 
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of some key cosmological parameters. Before reionization, most of the baryonic matter in 
the Universe is in the form of neutral hydrogen. The first gravitational bound objects are 
the collaps ed virialized da rk matter mini-halos where this neutral hydrogen resides. These 



mini-halos (Il77l : ll34l : ll35b provide a fluctuating background with a characteristic scale that 



can be used for weak lensing studies. 

For weak lensing studies, one would like to have a background at a high source redshift 
that also exhibits structures at small scales. Clearly the CMB satisfies the first criterium but 
it is smooth at small scales. It was pointed out that the Universe at the epoch of reionisation 
might be the most natural place to look for applications of weak lensing studies. Clearly 
not only it is at a very high redshift, but also it emits brightly in the hydrogen hyperfine 
transition line and has structures on many scales ranging from several arcminutes to under 
a milliarcseconds. Experiments which will target these particular redshifts and wavelengths 
are being planned, including some that are already undergoing construction, e.g. PAS 
LOFAR (1156F1 Other experiments which will have low signal-to-noise include T-RE 



3 



and CATWAL^. 



5 



Presenting a detailed calculation of signal to noise analysis (Si) argues that weak lens 



ing studies of epoch-of-reionization gas can constrain the projected matter power spectrum 
to very high accuracy. Use of tomography can further increase the level of accuracy with 
which certain cosmological parameters can be constrained. It was claimed that with such a 
technique, the neutrino mass could be constrained with an accuracy of 0. 1 meV. The infla- 
tionary gravity-wave background and consequently the inflationary dynamics as encoded 
by the Hubble parameter during inflation can be constrained with high accuracy too. These 
calculations show that such 21cmm weak lensing observations are an order of magnitude 
better than those from galaxy surveys. However such optimistic scenarios will require re- 
solving each of the 10^^ mini halos that will be observed on the sky. 

Clearly several problems ar ise w hen one tries to map gravitational lensing at such an am- 



bitious scale. For example (12071) have shown that synchrotron emission from ionised gas 



can outshine the 21cm radiation. However, pointed out that such components can be 
removed by power-law spectra from spatial fluctuations. 



^ http://astrophysics.phys.cmu.edu/ jbp/pas t6.pdf 

^ http://www.Iofar.org/ 

^ http://orion.physics.utoronto.ca/sasa/DownIoad/poster/casca_poster.pdf 

^ ftp ://ftp. astro. unm.edu/pub/users/john/AONov03 .ppt 
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9 Weak Lensing of the Cosmic Microwave Background 



Observation of cosmic mic r owav e ba ckgr ound (CMB) radiation is one of the cleanest 
probes of cosmology (|126l : 11271 : l63l : 11291 : H 549). However lensing of CMB photons 
by intervening mass clumps can provide additional information about the structure and dy- 
namics of the Universe. Besides, lensing of the CMB provides information at larger scales 
and higher redshift than can be reached by any other astronomical observations. Weak 
lensing of the CMB is responsible for many obse rvable effects which have been studied in 
extensive detail - for a detailed review see (Il68h . Calculations o f these ef fects have been 
made both for temperature and polarisation anisotropics (|31; 49; 172 : 254 : 193|) . 



9.1 Effect of weak lensing on the temperature and polarisation power-spectrum 



Lensing broadens the acoustic peaks and enhances power at small angular scales. These 
features are non degenerate with the standard cosmological parameters. This allows a de- 
termination of the lensing amplitude by comparing the observed spectra with CMB spectra 
of different mass fluctuations. The magnitude of distortion is sensitive to the level of mass 
fluctuations as a function of redshift and scale. This in turn depends on the background 
cosmology. In this way lensing of CMB can be used as a window to probe the fluctuations 
in the dark matter distribution over a huge range of redshifts and length scales. 

Lensing generates a non zero B-mode polarization from a purely E-mode. (Note that here 
we are referring to E and B modes in the CMB polarisation map, not in the shear map). 
B-modes generated during inflation by tensor mode perturbations or gravitational waves 
produce a distinct spectrum which can be used to discriminate between different classes 
of inflationary models. Lensing of E-modes produces B-modes which dominate the more 
primordial signal on small scales, and which are considered mainly as a source of confu- 
sion in polarisation experiments. Detection of a primordial B-mode signal in the presence 
of lensing therefore pushes observational strategies towards favouring larger sky coverage 



I55I : 12551) . A low sky coverage and presence of boundaries also causes additional con- 
fusion by introducing mixing of E and B modes. It is however useful to note that as the 
inflationary B-mode signal is mainly significant on large angular scales experiments with 
low resolution can also be very promising. High resolution B-mode detection experiments 
will typically employ a "delensing" step in data reduction to effectively restore the unlensed 
sky from the lensed data. 

Typically the power spectrum of the lensed map is computed by a series expansion in the 
deflection angle. A perturbative expansion in the deflection angle is a transparent way to 
understand most of the lensing effects. Non-perturbative evaluations of lensing effects are 



carried out by considering the correlation function (|254l : l45l : |5 ll : |3 1 ll : l55l : Il78j) . 
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Fig. 21. Left Panel: Power spectrum for the temperature anisotropies in the fiducial ACDM model 
with X = 0. 1. In the case of temperature, the curves show the local universe contributions to CMB 
due to gravity (ISW and lensing) and scattering (Doppler, SZ effects, patchy reionization). See, 
Cooray,Baumann & Sigurdson (2004) for a review on large scale structure contributions to temper- 
ature anisotropies (Figure Courtesy: Asantha Cooray). Right Panel: CMB B-mode polarization. The 
curve labeled 'IGWs' is the IGW contribution with a tensor-to-scalar ratio of 0. 1 with (solid line; 
T = 0.17) and without (dashed line) reionization. The curve labeled 'lensing' is the total lensing 
confusion to B -modes. Thin lines show the residual B-mode lensing contamination for removal of 
with lensing out to Zs- See, Sigurdson & Cooray (2005) for details. 

9.2 Non-Gaussianity in the CMB induced by Weak Lensing 



Primordial perturbations produced during inflation are very nearly Gaussian and the linear 
evolution of these perturbations up to the last-scattering surface does not generate any non- 
Gaussianity. Weak lensing by the interve ning matter distribution however does introduce 
non-Gaussianity in maps of the CMB sky (l31l|) . 



All odd-order correlation functions vanish for a Gaussian random field. On the other hand 
weak lensing generates a non-zero contribution through correlations between the large- 
scale temperature and lensing potential. This cross-correlation which is due to the In- 
tegrated Sachs-Wolfe (ISW) effect at large angular scales can therefore be a very use- 
ful tool to pr obe the growth of perturbations and the expansion history of the Universe 
(11241 : l256l : ISSi : l82l: l83l: 1841). This implies that such studies can be v ery u seful tools for study- 
ing the dark energy equation of state or the neutrino mass (11661) . These correlations 
are negligible on smaller scales. Besides non-linear evolution of the matter distribution 
as well as late time non-linear effects including the Sunyaev-Zeldovich (SZ) effect also 
contribute to the bispectrum. The four-point correlation function is the lowest order non- 
zero correlation function which does not vanish in the absence of any the cross -correlation 
between the low redshift mass distribution and the CMB temperature distribution. The non- 
Gaussianity studies involving polarisation fields are very similar to the temperature case, at 
least for the E-mode polarisation. Various combinations of four-point correlation func tions 
involving E- and B- mode polarisations have been studied in the literature (j23|; 11221) . For 
future high resolution all-sky polarisation surveys such studies will be feasible. 
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9.3 Weak tensing effects as compared to other secondary anisotropics 



The thermal Sunyaev-Zeldovich (tSZ) effect, which is caused by the scattering of photons 
from hot electrons in clusters, is a dominant anisotropy contribution on small scales. Due 
to a very characteristic frequency spectrum such a component can however be readily sepa- 
rated from primordial anisotropics. Another secondary contamination, the kinetic Sunyaev- 
Zeldovich (kSZ) effect, has the same frequency dependence as the primary anisotropics, so 
is harder to separate. Similarly, lensing does not cause any frequency shift in an otherwise 
perfect blackbody spectrum of primary CMB. Thus non-linear sources, such as kSZ are po- 
tential sources of confusion, when trying to understand the effect of lensing on the CMB. 
Current uncertainties in the reionization history and topology of reionization patches also 
make it more difficult to model. For polarisation spectra the kinetic SZ is sub-dominant, 
simplifying the situation to some extent (,289i:l80l:l43l:l262l). 



9.4 Lensing of the CMB by individual sources 



On small scales, the CMB lacks power. As mentioned earlier, most of the power is gener- 
ated by secondary anisotropics and transfer of power from larger scales to smaller scales 
due to lensing. This raises the possibility of detecting individual cluster-mass objects in 
CMB maps by their effect of lensing on the CMB. Unlike other probes such as the SZ 
which depends on baryon physics , suc h studies can constrain the physical mass distribu- 
tion of individual objects directly (| 1671 : 1641 : 1531). 



9.5 Future Surveys 



In the future, satellite experiments such as Plancl<[Z] (to be launched in 2008) will be in a 
good position to detect the effect of lensing in temperature power spectra.In the case of 
polarisation, experiments such as CLOVERL^I and QUIE10, which are bei ng planned, will 
have detection of B-mode from inflation as their primary science driver (|169l : l200l : l269l : 
171|) . Ongoing experiments such as may reach the required sensitivity, but for all 



of these, lensing is the main source of confusion on small scales for such experiments. 



^ http://www.rssd.esa.int/index.php ?project=Planck 

^ http ://w w w. astro. cf. ac . uk/groups/instrumentation/proj ects/clo ver/ 

^ http://quiet.uchicago.edu 

http://www.astro.cf.ac.uk/groups/instrumentation/projects/quad/ 
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10 Weak lensing and External data sets: Independent and Joint Analysis 

Weak lensing is a very powerful probe of the projected dark matter clustering, and in prin- 
ciple it is very good at constraining the cosmological parameters playing a dominant role 
in the structure growth. The main limitation comes from the strong degeneracy between the 
dark matter power spectrum normalisation and the matter density Q.jn (see the introduc- 
tion of Section [5] and Eq. (|2.9l) ). These two parameters are degenerate because any change 
in Q.^ can be balanced by an appropriate change in Og slowing down or accelerating the 
growth rate. In (l22t) it was shown that the shear variance at scale 0^ scales roughly as 

This equation, which is the first term of the perturbation series on the mass density contrast, 
assumes a power-law power spectrum with constant slope n and a single source redshift Zs- 
It has a very limited application but it has a pedagogical value in that it shows that the pa- 
rameter degeneracy between Og and Q.^ extends to the source redshift as well. Therefore, 
weak lensing in 2D can become a high precision cosmology tool only if it is combined 
with another cosmology probe and provided we have a good knowledge of the source red- 
shift distribution. The simultaneous observation of non-linear and linear scales in lensing 
surveys shows some features in the projected mass power spectrum which helps to lift this 
degeneracy (|136) . Alternative approaches which can also lift degeneracies were discussed 
in Section IH using 3D information, and in Section [51 using higher-order correlations be- 
yond the shear variance. However, like any other cosmology probe, this is not enough for 
weak lensing to be a high precision cosmology tool alone. In this Section we review the 
advantages of combining lensing with different probes and we outline the gain on cosmo- 
logical constraints in the future missions. 

10.1 With CMB, Supernovae and Baryon Acoustic Oscillations to probe Cosmology 

Weak lensing alone provides a measure of O^o^, provided that the redshift of the sources 
is known. This means that weak lensing best performance is in the measurement of the 
amplitude of the projected mass power spectrum. For this reason, it is very powerful at 
breaking the parameters degeneracies seen in other cosmology probes, which are usually 
sensitive to other combinations of cosmological parameters. It is known, for instance, that 
a precise measurement of (5% and Q.^ can be obtained from the combination of lensing 
and CMB (298,; .50[). Figure |22] sho ws the set of parameters to be combined for an optimal 
joint CMB-lensing analysis fisj- This result, obtained for the CFHTLS and WMAPl 
surveys, is not modified for a different choice of lensing and CMB data sets. Note that 
the best improvement is obtained on parameters which the shape of the dark matter power 
spectrum is the most sensitive to: the mass density Qm, the power spectrum normalisation 
Og, the reduced Hubble constant h, the primordial power spectrum slope n^ and the running 
spectral index a = — dln(n,s)/dlnA;. 
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Fig. 22. Figure showing the most orthogonal parameter degeneraties between CMB (WMAPl) and 
weak lensing (CFHTLS). The set of parameters is Q.^^, ag, cOc = 0.ch^, rig, and h (from Tereno et 
al. 2004). 



For the lensing signal alone, an improvement of a factor ~ 3 in the parameter errors is 
expected if the lensing signal is combined fo r diff erent source redshift slices instead of 
measured from the broad source distribution (|276l). This is the well known tomography 
technique which requires the measurement of photometric redshifts (see Section |431) . Un- 
fortunately, for these parameters, the improvement cannot be further increased with a larger 
number of source s lices : one can show that the signal-to-noise saturates for a maximum of 
3-5 redshift slices (|277|) . which is a consequence of the fact that the lensing selection win- 
dow is a rather flat function of redshift where the signal is the strongest. Anot her so urce 
of improvement for the lensing constraints is the use of higher-order statistics (( 276h . see 
Section [5] of this review): the probe of the non-linear regime of structure formation helps 
considerably in the determination of the precise moment when large scale structures be- 
come non-linear and at which scale. T his ev ent is a strong function of the cosmological 
parameters, in particular of d ark energy (Il83l : ll8[) . Some measurement of high order statis- 
tics have been done (|26l : l213h . but this area of research is still in its infancy. 
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Fig. 23. Panels showing the measurement forecast of the dark energy equation of state parameter Wp 
and dark energy density D-de from the baryon oscillations (top-left panel), supernovae la (top-right 
panel) and weak lensing (bottom-left) from the DETF report (2006). Wp is the equation of state 
parameter of the Dark Energy parameter at an intermediate 'pivot' redshift, usually a little less 
than 0.5. Solid and dashed lines are for optimistic and pessimistic surveys respectively. The surveys 
characteristics are described in the text. The bottom-right panel shows the constraints from the 
PLANCK CMB experiment with a lensing survey (Takada & Jain 2005). 

The main science driver for cosmology has become the measurement of the Dark Energy 
equation of state. The Dark Energy Task Force report (DETF) (2) and the ESA-ESO work- 
ing group report (jlJLOj) are summaries of where cosmology is heading to for the next decade: 
there is a consensus among cosmologists that the goal of measuring the dark energy param- 
eters can only be achieved from the joint analysis of several cosmological probes. Figure 
|23] shows the dark energy parameters forecast for the future experiments (which includes a 
pessimistic and optimistic cases). Most of the panels on this figure are extracted from the 
DETF report, and they show the various degeneracies for different cosmological probes. 
The three projects that DETF has considered are the following: 

• Baryon Oscillations: 20000 square degrees, ground-based survey. Photometric redshifts 
cover the 0.2-3.5 range, and their precision is 0.01 for the optimistic case and 0.05 for 
the pessimistic one. 

• Supernovae: ground-based survey with 300000 supernovae. The photometric redshift 
accuracy is 0.01(1 +z) for the optimistic case and 0.05(1 +z) for the pessimistic. 

• Weak Lensing: 20000 square degrees, ground-based survey. The shear calibration is 
fcal = 0-01 and photometric redshift accuracy — 0.01(1 +z) for the pessimistic case. 
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fcal = 0.001 and = 0.001(1 +z) for the optimistic case. 



Note that these surveys are all providing predictions of very accurate determination of the 
dark energy equ ation of state Wp and energy density Q.^^, far better than any joint analysis 



would do today (|267l) . Here Wp is the equation of state at a n inte rmediate reshift, typically 



around 0.4. The bottom right panel on Figure [23] is from (12761) and shows the joint con- 
straints on Wp and Q.^^ from the future CMB PLANCK mission and a lensing survey cov- 
ering 4000 square degrees, combining the two and three-points statistics. Figure [23] shows 
that the supemovae, CMB and lensing (or BAO) have very different degeneracies in the 
dark energy parameter space, which offers an optimal complementarity in the combination 
of cosmology probes. It is believed that only the combination of all probes together will 
provide convincing constraints regarding the dark energy, although each individual probe 
seems to predict very accurate results. The main limitation being the systematics, it is in- 
deed important to have more constraints than parameters we wish to measure: weak lensing 
for instance may appear as the most powerful probe. However, the lensing signal is rather 
featureless compared to CMB and BAO, which makes it more vulnerable to systematics 
like shear calibration and photometric redshift inaccuracies (see Section [6]). The future of 
cosmology certainly lies in a joint analysis of surveys, and lensing surveys play a particular 
role in the sense that this is the only probe which can provide an unbiased measurement 
of the dark matter fluctuations amplitude. Note that the constraints on quintessence mod- 
els can be improved by the combination of various cosmology probes, such as lensing and 



SNela as demonstrated in ( 230 ). 



The success of the joint analysis is subject to the validity of our Cold Dark Matter and 
Dark Energy picture. Given the unknown nature of these ingredients, it is possible that our 
description might be incomplete or wrong. We must therefore keep an open mind and con- 



sider alternative interpretations of the data. The modified gravity theory proposed by (11961) 



lacks a solid physical moti vatio n and it appears to be unable to explain the halo flattening 



seen in weak lensing data (1115[) . Nevertheless, this study initiated the idea that a modified 
theory of General Relativity could solve the dark matter and dark ene rgy p roblems. These 
models include the string-motivated braneworld scenarios (69). Ref.(286) suggested that 
modified gravity can be tested directly from a verification of the Poisson equation 

= AnGp. (10.2) 

The gravitational potential 4> can be measured from weak lensing by large-scale structures, 
and the mass density p from the galaxy distribution, which at large enough scale can be 
assumed to be an unbiased tracer of the mass distribution, as shown in (|30lb from the 2dF 
Galaxy Redshift Survey. 

A general problem with modified theories of gravity is the significantly inc rease d level of 



degeneracy due to a larger number of degrees of freedom. As pointed out by ([131|) . it is nec- 
essary to measure independently the growth rate of structures and the Universe expansion 
rate in order to reduce the degeneracy to a reasonable level. Modified gravity also might 
not be the cause of all the dark components of the Universe. The investigation of a realistic 



mix of dark component and modified gravity in (|131|) shows that neglecting the modified 
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gravity severely biases the dark energy parameters, even with a joint Supemovae, CMB 
and weak lensing constraints. To make things worse, the addition of inflation parameters 
introduces quantum correction, the running spectral index, to the primordial mass power 
spectrum which is a scale-dependent effect. Yet another degree of complication comes from 
the small-scale amplitude prediction which is highly non-linear. Dark Matter itself could 



also be self- interacting as suggested by (|40|) . resulting in density-dependent observational 



effects (I317h . 



The present situation is that the simple Big-Bang scenario (i.e. no alternative theories) will 
indeed be accurately constrained by a joint analysis of different surveys. A general scenario 
including alternative theories might be harder to constrain. Whether or not this is possible 
with the next generation of surveys remains to be demonstrated, but it is clear that weak 
lensing by large scale structure plays a central role in being sensitive to the "dark matter" 
whether it is real dark matter or modified gravity. 



10.2 With Galaxy Surveys to probe bias 



Weak lensing is the only reliable technique able to probe the dark matter distribution up 
to redshift of a few using optical surveys, or at even higher redshifts from lensing of the 
CMB and the 21cm line. The peculiar velocity field, which is also an unbiased tracer of 
the matter distribution, is no longer accurate at distances larger than a few hundred Mega- 
parsecs. The combination of lensing with the galaxy distribution is therefore a unique way 
of constraining the relative amount of matter with respect to light, the so-called bias. This 
can be done out to reasonably high redshift, giving astronomers access to the dark side of 
galaxy formation over a wide range of its evolutionary history. 



10.2.1 Galaxy biasing 

An apparent limitation of lensing surveys is that the mass is only seen in projection, imped- 
ing a 3-d imen siona l probe of the bias. A method for alleviating this problem was proposed 
by ('292') and (^236^. It is an alternative to the 3D mass reconstruction technique discussed 
in Section m which was designed specifically for potential or density measurements. We 
assume a population of foreground galaxies with a known narrow redshift distribution nf{z) 
centered on Zf (in units of arcmin^), from which we want to measure the mass-to-light ra- 
tio. The lensing signal is measured from a background galaxy population which can have 
a broad redshift distribution nb(z) (normalized to unity), while the technique discussed in 
Section m requires an estimate for all redshifts. The cross -correlation of the lensing signal 
with the foreground galaxy density distribution provides a measurement of the bias at red- 
shift Zf. The scale dependence of the bias can be obtained by filtering the lensing signal and 
the galaxy density with an aperture filter, which is a narrow band filter (Fig.[2l). 

Assuming that the galaxy number density contrast 5gai is proportional to the matter density 
contrast 5 with the proportionality factor defined as the bias parameter b, the lensing-galaxy 
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density cross -correlation at scale 0j is given by (|292|) : 



see Eqs.(l2:9l)- (l2l2l ).(l34l) and (l3TT0l) in Sections [21 and [3l Here the weight w{x) along 
the line of sight depends on the redshift distribution n\j{z) of the background sources as 
in Eq. (|2.10l) whereas hf{x) is the redshift distribution of the foreground galaxies (with 
^f(x)dx = nf{z)dz) of "number counts" 9{_ with respect to the aperture of radius 0^,. P{k;x) 
is the time-evolving 3-D matter power spectrum (where % is a parameterization of the red- 
shift along the line of sight) and W^^^{iBs) is the square of the Fourier transform of the 
aperture filter, which peaks at some effective wavelength £eff ~ 5/05, see Fig. [2l For a 
narrow foreground redshift distribution, the function hf{x) peaks at some radial comoving 
distance Xf{zf) ■ Eq. (|10.3l) shows that the cross-correlation is dominated by P (/ceff ; Xf ) where 
^eff ~ 5/(© (Xf)0i). In order to extract the bias parameter b, we need to define the variance 
of the number density fluctuations for the foreground galaxies: 

W^(e„)) = ft^(^e?)^/ dx|^ l^iP (10.4) 

and define the ratio R 

„_ (M.p(e,)at(ec)) 



It was shown in (|292f) that this ratio is nearly independent of scale for all cosmologies 
and any dark matter power spectrum, unless the bias parameter is scale-dependent (this 
statement was shown to be valid in the non-linear regime as well). The bias at angular scale 
k^fl and redshift Zi can therefore be measured from weak lensing data. 



Refs. (lll3l) (12641) have performed the only application of this technique. (11131) measured 
the bias from a combination of the VIRMOS (298) and RCS (111) surveys, which is shown 
on Figure [24l This analysis shows a significant scale dependence of the bias, although its 
calibration is still uncertain due to incomplete knowledge of the source redshift distribution. 
With the GaBoDS surveys, (|264t) found = 0.8 ±0. 1 and r = 0.6 ±0.2, which is consistent 
with dlBl) . 



10.2.2 Galaxy-galaxy Lensing 

The probe of galaxy biasing can be extended dow n to g alactic halo scales with a technique 
called galaxy-galaxy lensing. First proposed by (|285b and then formalised by (35), the 
idea is to cross-correlate the shape of distant lensed galaxies with foreground galaxies. 
As shown in the left panel of Figure [251 the background galaxies lensed by a foreground 
galactic halo are preferentially tangentially aligned with respect to the foreground galaxy. 
The situation is identical to lensing by a cluster of galaxies, but the lensing by individual 
galaxies is much weaker, since the amplitude of the effect scales as the mass of the lens. 
To measure the galaxy-galaxy lensing, one must therefore stack the lensing signal behind a 
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Fig. 24. Galaxy biasing constraints from the VIRMOS and RCS weak lensing surveys (Hoekstra et 
al. 2002). The left panel shows the cross-correlation coefficient r and biasing b as function of scale. 
The right panel, top plot, shows the ratio b/r and its prediction for a concordance model = 0.3, 
Ha = 0.7. The bottom plot shows the most convincing test of residual systematics obtained by 
rotating the lensed galaxies by 45 degrees, which is supposed to cancel the lensing signal. 

large number of foreground lenses. The average shear as function of angular distance from 
the center gives an estimate of the average halo profile around the foreground galaxies (see 
right panel in Figure [25]). 

The tangential shear of a lensed galaxy at position angle with respect to the lens on the 
sky is given by 

Yf = -ei cos(20) + 62 sin(20). (10.6) 
It is straightforward to show that, when averaged over all position angles, the mean tan- 
gential ellipticity is unchanged if a constant is added to the galaxy ellipticity e = (ei,e2)- 
For this reason, galaxy-galaxy lensing is robust against an imperfect Point Spread Function 
(PSF) anisotropy correction, as long as the latter is a slowly-varying function of position. 
The lensing signal is therefore relatively easy to measure, even if the amplitude of the signal 
is low. 

Given that the foreground galaxies span a large range of velocity dispersion and luminosity, 
scaling relations are needed in order to calibrate the expected lensing signal to the same 
fiducial galaxy with luminosity L^, and size s^,. If o^, is the velocity dispersion of the fiducial 
galaxy, the scaling relations are: 



L 



o 



(10.7) 



where the latter corresponds to the TuUy-Fisher or Faber- Jackson relation for spiral and 
elliptical galaxies respectively if T] = 1/4. (L, 5, a) are the luminosity, scale and velocity 
dispersion of one foreground galaxy. The lens mass model chosen by (jishjs a truncated 
isothermal sphere (TIS), which is also frequently used by several authors (|l30l : lll5|) . The 
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Fig. 25. Left panel: schematic illustration of galaxy-galaxy lensing. Right panel: Galaxy-galaxy 
lensing signal from the GEMS survey (Heymans et al. 2006). The soUd line shows the best fit NFW 
mass profile, assuming an average lens redshift of zi = 0.65. 



mass density of the TIS is given by 



P(r) 



(10.8) 



whose total mass Mtot is (|115|) 



Mtot = 1.3 X lO^^h^^ Mq 



lOOkms 



-1 



IMpc 



(10.9) 



The tangential shear is calculated from the mass model (e.g. Eq llO.81) . and can be compared 
to the data provided that an estimate of the lens and source redshifts is know. In general, the 
free parameters measured with galaxy-galaxy lensing are the fiducial velocity dispersion o^, 
and truncation radius of an galaxy. The apparent magnitude of each foreground lens is 
used to estimate its absolute luminosity L (which may require an appropriate ^-correction 
to be included), then the velocity dispersion a is o btained from the scaling relations. The 
maximum likelihood technique developed by (l235b ensures an optimal analysis, in particu- 
lar if individual photometric redshifts can be obtained. Note that simple mass models such 
as the singular isothermal sphere (SIS) do not have a truncation radius: 



(10.10) 



in which case one could constrain the scaling re lation s, Eq. (ll0.7l) . like the parameter r\, 
in addition to the fiducial velocity dispersion o^,. ( 130l) and (1 162) measured a r\ parameter 
very close to the TuUy-Fisher and Faber- Jackson relation. 



At low redshift, the most extensive galaxy-galaxy lensing analysis was performed by (19 11) 
on the SDSS data set. They found a Virial mass M200 = 5 - 10 x 10^^ H^^Mq for = 
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10^^ h ^Lq, depending on the galaxy color and morphological type. In the redshift range 
0.2 — 0.7 a Virial mas s of M200 — 4 — 8 x 10^^ h~^MQ, the less massive corresponding 



to the bluest galaxies (|101l : 11621) . (|259|) measured the galaxy-mass bias from the SDSS 



galaxy-galaxy lensing signal, and found a constant bias over correlation coefficient ratio 
Z?/r=1.3±0.3for = 0.27, which is in agreement with the cosmic shear study presented 
in Section [ro.2. II (see Figure l24)). 



Recent studies use the Navarro-Frenk- White density profile as the parametric mass model, 
which is particularly well suited for galaxy-galaxy lensing in mo re massive structures such 
as galaxy groups (|ll2h and galaxy clusters (12051 : |79|) . In (|l8lb the authors measured the 
lensing around Luminous Red Galaxies (LRG) from the Sloan Digital Sky Survey. The 
LRG sample was split in a bright and faint sample at Mcut = —22.3. LRG are particularly 
good foreground targets for probing larger mass halos because they are known to be present 
in the core of groups and cluster of galaxies. Figure [26] shows the measured signal against 
several mass models. It is particularly interesting to note that the flattening of the dark 
halo profile is clearly visible, leading to a con centr ation parameter of c ~ 5 — 7, in perfect 
agreement with Cold Dark Matter predictions (|206h . 



(Il01|) marginally measured a halo over lens ellipticity ratio of e/,/eg ~ 0.8 ± 0.2, consistent 
with the Cold Dark Matter Scenario (|68[) . A significant, but contradictory, measurement is 
shown in(^, who found eh/ eg = 0.1 ±0.06 for red galaxies and 0.8 ±0.4 for blue galax- 
ies. The baryon fraction can also be measure d by galaxy-galaxy lensing from a comparison 
of the lensing to the stellar mass. (U6) and (|l04|) find a virial to baryon mass-to-light ratio 
between 50 and 100, and they show that massive early type galaxies have a stellar to total 
baryon fraction of ~ 10%, which indicates that these galaxies are not efficient at producing 
stars. 



10.3 With Sunyaev-Zeldovich studies to probe small scale baryonic physics 



The energy of the CMB photons is boosted by scattering on the free electrons contained in 
the plasma of the Intra Cluster Medium (ICM). It is the source of the Sunyaev-Zeldovich 
(SZ) effect, which shifts the spectral energy distribution of CMB photons for the lines-of- 
sight containing hot, ionized, gas. The SZ effect is a probe of the hot baryon distribution, 
which can be compared to stellar and lensing mass distributions in order to learn about the 
cluster physics. If one assume that the gas is isothermal, SZ is a measure of the electronic 
density ne(9,z) projected along the line of sight: 

-> kT C - 
j(e) = ^C5t dzn,id,z), (10.11) 

where Oj- is the Thompson scattering cross-section and T the temperature of the gas. The 
SZ effect is therefore independent on the cluster redshift, as opposed to the lensing effect 
which has its maximum sensitivity at mid distance between the observer and the sources. 
Lensing and SZ are complementary because they are both linear in the density: for this 
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Fig. 26. Galaxy-galaxy lensing signal for the Luminous Red Galaxies (LRG), used as tracers of 
galaxy groups and small clusters (Mandelbaun et al. 2006). 43335 LRG from the Sloan Digital 
Sky Survey were used. The left and right panels correspond to faint (Mr > —22.3) and bright 
{Mr < —22.3) samples respectively. 

reason they probe the same regions of galaxy clusters, while X-rays for instance scale as 
the density squared, which is more sensitive to the cluster core. 



One approach of the combined SZ and lensing data sets consists of working on individual 
clusters. Both data sets can be used to predict the X-ray emission of the cluster which 
provides a direct test of the cluster dynamical equilibrium. The 3-dimensional halo shape 
(assuming axial syn imet r y) ca n be re constructed from the combination of SZ, lensing and 
X-ray observations 13H217|) . (|25 li) have shown that large lensing and SZ surveys (a few 
hundred square degrees) would be able to detect the evolution of the mass-SZ luminosity 
relation with redshift, which in turn is a measure of the cluster baryonic physics. This is 
particularly important for the understanding of the source of energy maintaining the tem- 
perature of the intra-cluster plasma at high temperature. The other approach is statistical: 
the goal is to use the cluster number counts to probe cosmology (16), the lensing data is 
used to measure the cluster masses, but this requires a full sky survey with a mass sensitivity 
down to a few 10^"^ Mq to provide interesting constraints. 



One should remember that cluster masses derived from lensing are subject to large noise 
and projection effects (IllOl : although recent studies seem to s how t hat the discrep- 
ancy between dynamical, lensing and X-ray masses is rather small The ideas de- 
veloped around the combination of SZ and lensing yet remain to be put into practice: in- 
struments such as the Cosmic Background Interferometer, Atacama Cosmology Telescope, 
Arcminute Imager and South Pole Telescope, all with angular resolution of the order of a 
few arcminutes, are very promising for this purpose. The understanding of cluster physics 
from SZ and lensing is not without any consequences for cosmological studies: statistical 
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lensing for instance probes the projected mass power spectrum down to arbitrarily small 
angular scale (if the statistical noise is low enough). At the cluster scale, typically one 
arcminute, we know that baryon cooling and heating modify the c luster gravitational po- 
tential well, and therefore the dark matter distribution itself. ( 314 ) and QOS) have shown 
how the cluster physics could affect the power spectrum by ~ 10% below one arcminute. 
Taking into account this effect might be necessary for the next generation of high precision 
weak lensing studies. 



10.4 Weak lensing of supernovae and effects on parameter estimation 



Type la supernovae (SNela), which are believed to be the thermonuclear explosion of an 
accreting white dwarf, are standard candles with a small intrinsic dispersion around their 
average luminosity. Moreover, this dispersion can be further reduced to about 0.12 mag 
using an empirical relation between the peak magnitude and the width of the light curve 



(l223i) . This makes SNela excellent tools for observational cosmology. In particular, by 



measuring their apparent magnitude we can derive their distance from us and obtain the 
redshift-distance relation up to z ~ 1 which provides useful constraints on cosmological 
parameters (issb . This method has su pplied the main contribution to the discovery of the 



present acceleration of the Universe (12241 : |215|) and it is the basis of future proposals to 



probe the nature of dark energy through its equation of state (e.g. the SuperNova Accel- 
eration Probe, (0)). On the other hand SNela, like all radiation sources, are affected by 
gravitational lensing effects which can magnify or demagnify their observed luminosity. 
This is a source of noise for studies which intend to measure the redshift-distance relation 
but this effect may also be used by itself to constrain cosmology in the same manner as 
gravitational lensing distortion of distant galaxies. 

There are several important differences between SNela and galaxies as probes of gravi- 
tational lensing effects. Firstly, since SNela are point sources one uses the magnification 
(associated with the convergence k) rather than the shear y which we focussed on in previ- 
ous Sections. In particular, the magnification can be written as: 

t^=T-. i — whence ;U~1+2k for |k|<1,|y|<1. (10.12) 

(l-K)2-|y|2 

Second, whereas weak gravitational lensing only modifies the observed ellipticities of 
galaxies at Zs = 1 by less than 10%, so that one needs many galaxies to extract the signal, 
the magnification of a type la supernova at Zs = 1 by gravitational lensing is of the same or- 
der as the intrinsic magnitude dispersion. Hence the signal-to-noise ratio is larger for SNela 
but since we have many more observed galaxies than SNela, accurate weak gravitational 
lensing effects have only been measured from galaxy ellipticities so far. Third, in order to 
derive the coherent shear on large scales one must cross-correlate the observed ellipticities 
of many distant galaxies over a window radius of the order of a few arcmin (since galax- 
ies are not exactly spherical). This leads to observables such as the aperture-mass Map or the 
mean shear over scale 0^ which probe matter density fluctuations over scales of the order of 
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Fig. 27. Left panel: the probability distribution 'P {/u) of the magnification within a ACDM uni- 
verse for SNela at redshifts Zs = 0.5 or = 1. From Valageas (2000a). Right panel: the constraints 
on cosmological parameters obtained from the dispersion of SNela magnitudes from a SNAP-like 
survey. The two sets of contours correspond to analyzing the data assuming the true non-Gaussian 
distribution (shaded) or a Gaussian distribution (unshaded). From Dodelson & Vallinotto (2006). 



'DQs- On the contrary, the magnification of each SNela only probes the density fluctuations 
along its line-of-sight, which corresponds to no smoothing (9^ = 0). Then, by measuring the 
probability distribution of observed SNela magnitudes (or its variance) rather than cross- 
correlating different SNela, one can probe the statistics of the convergence K with 0^ = 0. 
This means that the signal is dominated by scales where £^/k(^) is maximum, which are 
set by the matter power- spectrum. For CDM power-spectra this corresponds for sources at 
Zs = 1 to wavenumbers k ~ lOh Mpc^^ whereas smoothing galaxy ellipticites over 1 arcmin 
mainly probes smaller wavenumbers ~ 1/z Mpc^^ (287). Thus, weak lensing magnifica- 
tion of SNela allows us to probe density fluctuations on smaller scales than with galaxy 
ellipticities. This also implies that non-Gaussianities are more important for SNela grav- 
itational lensing distortions. On the other hand, this difference between the scales probed 
by galaxy shear maps and SNela magnitude distortions means that both effects are only 
weakly correlated so that weak-lensing shear maps obtained from surrounding galaxies are 



not very efficient to correct the SNela luminosities (|62|) . 



The rms magnification of SNela was computed by analytical means in (|76|) who found 
that this would not significantly decrease the accuracy of the distance-redshift relation at 
Zs < 0.5 used to measure the current acceleration of the universe but it could have a sig- 
nificant ef fect a t higher redshifts Zs> 1 , in agreement with the numerical simulations per- 



formed in (|303[) . In particular, although future surveys covering more than a few square de- 
grees will be unaffected, pencil beams surveys (< 1 deg^) suffer significan t con tamination 



(|57|) . Turning to the use of SNela as a tool to detect gravitational lensing, (W) computed 
analytically the rms magnification of SNela to find that in a ACDM universe one needs at 
least 2400 SNela at Zs = 0.5 (or 1 10 SNela at Zs = 1) to detect weak lensing from their ob- 
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served ma gnitu des. The PDF of the magnification T {jS) was computed from a hierarchical 
model in (287') as well as its impact on the measure of cosmological parameters through 
the distance-redshift relation. In particular, this study shows the strong non-Gaussianity of 
the magnification with an extended high-/j tail which follows the high-density tail of the 
underlying matter density field, as can be seen from left panel of Fig. [27l Then, the high- 
luminosity tail of observed SNela could be used to detect weak lensing since it should be 
significantly enhanced by gravitational lensing. However, the number of observed SNela 
is still too small to dra w de finite conclusions about a possible detection of weak lensing 



effects by this method (|306|) . 



On the other hand, since the amplitude of the gravitational lensing contribution to the dis- 
persion of observed SNela magnitudes depends on cosmological parameters it could be 
used to constrain cosmology (57). Thus, (65) found that 2000 SNela in the redshift range 
0.5 < < 1.7 should be able to constrain up to 5% the amplitude of the matter power- 
spectrum. However, they point out that one needs to take into account the non-Gaussianity 
of the weak-lensing magnification distribution in order to obtain correct estimates, as seen 
in right panel of Fig. [27] 

In order to eliminate the systematic errors associated with uncertainties on the intrin- 
sic SNela luminosity distribution and its possible dependence on redshift, it is possible 
to cr oss-correlate SNela observed magnitudes with foreground galaxies, as advocated in 



(|195l) . Indeed, in the absence of gravitational lensing this cross-correlation would vanish. 
Using a halo model (195) found that a gravitational lensing signal should be detected with 
~ 250 SNela at Zs = 1- This will be within the reach of future experiments (e.g. the SNAP 
satellite should observe thousands of SNela up to z ~ 1 .7) but current survey s are too small 



to detect a correlation between SNela magnitudes and galaxy overdensities (Il92|) 
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11 Summary and outlook 



Ten years ago, the detection of weak lensing by large scale structures was only a dream. The 
progress accomplished, summarized in this review, are remarkable. The main reason for 
this progress in essentially the development of the Charge Coupled Device detectors which 
had a major impact on the accuracy of galaxy shape measurement. Adequate detectors are 
probably more important than larger telescopes; weak lensing needs wide field of view and 
well sampled Point Spread Function. One sh ould remember for instance that the very first 
tentative to measure distorted galaxy shapes (|285|) failed because of inadequate technology 
detectors. 



In this review we mostly focussed on the statistical lensing, what one can learn on cos- 
mology from the shear or convergence field, rather then focussing on individual lenses like 
clusters of galaxies. Like most of the cosmology probes, the information is mainly encoded 
in the power spectrum. We have shown how different cosmological parameters affect the 
projected mass (convergence or shear) and how it can be measured. In particular a combina- 
tion of Og and Q.^ appears well constrained from current lensing surveys. Statistical lensing 
can also be used to probe the biasing between dark matter and light, which provides clues 
on the assembly of galaxies inside their hosting halos. Future surveys should be able to 
provide accurate measurement of the biasing history as function of scale. A particularly in- 
teresting aspect of lensing is that it can probe the non-linear regime without any assumption 
on how light traces mass, and therefore go beyond the traditional power spectrum analysis. 
This property gives access to a different sensitivity to the cosmological parameters, an im- 
portant feature to help breaking the parameter degeneracy between Q.^^ and Og, but it also 
probes the history of the gravitational collapse. In that respect, high-order lensing statistics 
can be used to test the role of gravity during the collapse. The tomography technique, which 
has just started to be applied to lensing data, is very promising in probing the Dark Energy 
equation of state. The redshift slicing of the lensed sources, or 3D analysis of the shear 
field, are the only ways to measure the growth rate of structures (from the power spectrum 
and bi-spectrum) which is very sensitive to the dark energy content of the Universe. An 
interesting alternative to this is the 3-dimensional reconstruction of the mass distribution 
which in addition will give us the distribution of the dark matter in space. The combination 
of statistical lensing with other cosmology surveys in other wavelengths was also shown 
to be very important for three reasons i) the sensitivity to the cosmological parameters is 
different and sometime orthogonal ii) the systematics can be drastically reduced (e.g. 21 
cm observation which offers perfect source redshift measurement) iii) this is the only way 
to probe the physics of the lenses beyond the simple mass-light relation. 

It appears that ray-tracing techniques are an essential ingredient of any precision weak 
lensing study. In particular this would be the only way to address many of the complicated 
higher-order lensing effects such as multiplanes deflections, source clustering, intrinsic and 
intrinsic- shear alignment which cannot be modeled with high precision. This is particularly 
relevant for future lensing surveys. Moreover, a proper assessment of the cosmic variance 
associated with the survey geometry can also be done from numerical simulation, which 
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is also important for non-linear scales (typically less than half a degree) where prediction 
from semi-analytical models is challenging. 



In 2006, weak lensing by large structures has just began to reach a status of scientific 
maturity: the first successes have demonstrated the feasibility of the technique and now 
begins an era of thorough weak lensing studies. The situation is similar to the Cosmic 
Microwave Background research in the pre-COBE era, before 1992: many independent 
groups have now measured the lensing fluctuation amplitude, and the next generation of 
lensing surveys will provide full sky, or nearly full sky, coverage. Table [5] summarizes 
the ongoing and future lensing surveys. The largest ongoing effort is the Canada France 
Hawaii Telescope Legacy Survey (CFHTLS). This survey represents a transition because 
this is the last one which requires a percent accuracy of galaxy shape measurement in 
order to be fully scientifically exploited. One percent accuracy is what we are currently 
capable of. All future surveys require a sub-percent precision level, which is still beyond 
our capability, as demonstrated by STEP. In fact, below the percent accuracy, many other 
effects will complicate the lensing measurement and analysis. 

All recent studies show that the only way to quantify these effects, such as intrinsic align- 
ment, shear-intrinsic alignment and source clustering, is to have a redshift estimate of each 
galaxy. This is only doable with photometric redshifts, which poses additional challenges: 
assuming one can get enough colors to obtain accurate photometric redshifts, there is no 
spectroscopic survey to help calibrating objects fainter than / ~ 24. The magnitude calibra- 
tion of brighter sources could be problematic, especially from the ground due to zero-point 
fluctuations for different wavelengths. An absolute zero-point variation is not a problem, 
but an explicit dependence with color could jeopardize photometric redshift estimates. It 
appears that the only way to obtain unbiased photometric redshift is to have uniform mag- 
nitude calibration and to conduct a deep spectroscopic survey in order to validate the tech- 
nique at faint magnitude/high redshift. Both might only be doable from space where pho- 
tometry is stable. The tunable laser project (HI) is an interesting solution for the zero-point 



issue, while the satellite GAIA (18 1|) . successor of Hypparcos, could help with absolute as- 
trometry calibration. Interestingly, it seems that the shape measurement problem is likely to 
be solved in the next two or three years, thanks to the Shear TEsting Program!^. The main 
limitation of weak lensing by large scale structures might therefore not be shape measure- 
ment anymore but multicolor photometric calibration, a very old astronomical problem! 



http://www.physics.ubc.ca/ heymans/step.html 
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Table 5 

List of forthcoming lensing surveys (adapted from Peacock et al. 2006). Surveys are sorted in three 
groups separated by a line. The top group show the current lensing surveys, the middle group shows 
survey starting in one year at the latest, and the group at the bottom is essentially not funded or 
partially funded projects. 



Survey 


Telescope 


Sky coverage 


Filters 


depth 


Deep Lens Survey 


CTIO 


7x4 deg^ 


BVRz' 


R=25 


CFHTLS-Wide 


CFHT 


170 deg^ 


ugriz 


jAB=24.5 


RCS2 


CFHT 


1000 deg^ 


grz 


jAs=22.5 


KIDS 


VST 


1500 deg^ 


ugriz 


jAB=22.9 


Pan-STARRS 


PSl 


30000 deg^ 


grizy 


jAB=24 


VIKING 


VISTA 


1500 deg^ 


zYJHK 


/ab=22.9 


Dark Energy Survey 


CTIO 


5000 deg^ 


griz 


/ab=24.5 


DarkCam 


VISTA 


10000 deg^ 


ugriz 


/ab=24 


HyperCam 


SUBARU 


3500 deg^ 


TBD 


TBD 


SNAP 


Space 


300/2000 deg^ 


Narrowband (0.35-1.6) 


TBD 


LSST 


6m ground 


20000 deg^ 


Narrowband (0.35-1.2) 


jAB=27 


DUNE 


Space 


20000 deg^ 


TBD 


jAB=25.5 



83 



Acknowledgements 



We thank Stephane Colombi, Asantha Cooray, Scott Dodelson, Wayne Hu, Bhuvnesh Jain, 

Rachel Mandelbaum, Yannick MeUier, Ue Li-Pen, Elisabetta Semboloni, Masahiro Takada 
and Alberto Vallinotto for use of their figures. It is a pleasure to thank Catherine Heymans, 
Sanaz Vafaei and Jonathan Benjamin for their help with some of the figures and tables. DM 
would like to thank Martin Kilbinger, Jerry Ostriker, Lindsay King, Tony Tyson and mem- 
bers of Cambridge Planck Analyis Center for many useful discussions. LVW is supported 
by NSERC, CIAR and CFI. DM is supported by PPARC. 



84 



A Analytical modeling of gravitational clustering and weak-lensing statistics 



In order to derive the properties of weak lensing observables, like the shear y, one first needs 
to specify the properties of the underlying density field. We briefly describe in this appendix 
two such models which can be used to predict weak lensing statistics, the hierarchical 
models presented for instance in (i29O0 and the halo model described in detail in (|52|) . 



A. 1 From density to weak-lensing many-body correlations 



At lowest-order weak-lensing observables can be written as linear functionals of the den- 
sity field, as seen in Eqs. (|2.9l) - (|2.16l) . Therefore, their many-body connected correlation 
functions can be directly writ ten in terms of the connected correlations of the 3D matter 
density field defined by (|212|) : 



^p(xi,...,Xp;z) = (5(xi,z)...6(xp,z))c 



(A.l) 



Note that for a Gaussian field we have = for p > 3. For a weak-lensing observable X 
defined as in Eq. (|2.11|) this gives in real space (|287|) : 



fXs /"^ p _ r y _ _, / 

(x^), = / dx nd5c, / Hde, \p ' 

^0 J —oo ■ J -1 



©01 ©02 



©0, 



(A.2) 



and in Fourier space (|290|) 



(Z^), = (27r)-2/-i UiwP [ f[d^j Wx{k^j^Qs) m±i) . ..Kk±p))c. (A.3) 



Here we used Limber's approximation ~ fc^ as for Eq. (l3.4l) and the longitudinal Dirac 
factor 5£)(fc||i + . . . + k\\p) has been factorized out of the correlation (5..5)c. Next, from 
the correlation functions one can obtain the full probability distribution function (PDF). 
Indeed, if we define the generating function cpx(j) of the cumulants {XP)c as: 



^x{y) = £ / with 



p=2 



(A.4) 



where we used (X) = 0, then one can show that the PDF Tx{X) is given by the inverse 
Laplace transform: 



(A.5) 



Thus, in order to compute (Px{X) one only needs to derive cp;c(j), which may be written in 
terms of the density field from Eq. (IA.2|) or Eq. (IA.3l) . 
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For the smoothed convergence ic the cumulants (k^)c correspond to averages of the density 
cumulants over cylindrical cells along the line of sight and they can be obtained with a 
good accuracy from the 3D density cumula nts average d over spherical cells. If we also use 



a mean-redshift approximation one obtains (|287l : l288l : ll3i) : 







|Kmin|^95 ( 


V 1 i^min 1 / 



with K^n = - J dxw, (A.6) 

where we introduced the minimum value K^in of the convergence, which corresponds to an 
empty line of sight (5 = —1, see Eg. (12. 101 )). This gives: 

I i^min I \ I ^-min | | ^-min | / 

where (p§ and !P§ are the generating function and the PDF of the 3D matter density contrast 
at the mean redshift and scale probed by the smoothed convergence ic (they depend on the 
angular radius 0^ and the galaxy distribution n{zs))- Thus, within this simple approximation 
the PDF of the projected density field ic is directly expressed in terms of the PDF of the 
underlying 3D density contrast. 



A.2 Hierarchical models 



For more intricate observables like the shear or the aperture-mass which involve compen- 
sated filters one cannot perform approximations such as Eqs. (IA.6l )- (IA.7l) and it is not possi- 
ble to approximate high-order cumulants or the PDF of weak-lensing observables in terms 
of those of the smoothed density contrast. Therefore, one needs to specify the detailed 
angular behavior of the many -body correlation functions ^p(xi , . . . , x^) . A simple prescrip- 
tion is_grovided by the general class of "tree-models" defined by the hierarchical property 



(12261: 189h 



^p(xi, . . . ,x^) = £2},'^^ £ n ^2(x„x,) (A. 8) 

(a) ta p— 1 

where (a) is a particular tree-topology connecting the p points without making any loop, 

(a) 

2p is a parameter associated with the order of the correlations and the topology involved, 
ta is a particular labeling of the topology, (a), and the product is made over the (p — 1) 
links between the p points with two-body correlation functions. Then, as seen in (128 8b 
the 2D correlations (Op involved in weak-lensing cumulants such as (IA.2I) exhibit the same 
tree-structure, with: 



(Op{Qu...,Qp;z)= rfldx^^pi^ _ ;z . (A.9) 

J-^i=2 \ ©01 'DQr ' 



(a) 

In order to perform numerical computations we need to specify the weights Qp ' . Thus, the 

(a) 

"minimal tree-model" corresponds to the specific case where the weights Qp are given by 
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(21): 



n 



(A. 10) 



vertices of (a) 

where Vg is a constant weight associated to a vertex of the tree topology with q outgoing 
lines. The advantage of this minimal tree-model is that it is well-suited to the computation 
of the cumulant generating functions as defined in Eq. (|A.4| ). Indeed, for an arbitrary real- 
space filter, F(x), which defines the random variable s as: 



|dxF(x)5(x) and ^ = (5^), 



(A.ll) 



it is poss ible to obtain a simple implicit expression for the generating function, ^s{y), see 
SiSil): 



x(x) = -3;|dx'F(xO 



x(x)C;[x(x)] 



(A.12) 
(A.13) 



where the function ^v(x) is defined as the generating function for the coefficients Vp: 



^v(x) = E —r 
p=i p- 



Vp with vi = 1. 



(A. 14) 



Since the 2D correlations (Op obey the same minimal tree-model we can perform the re- 
summation (|A.12| )- (|A.13I ) which yields (1251: IbI): 



(A. 15) 



where we introduced the 2D generating function cpcyi. associated with the 2D correlations 
(Op, given by the resummation: 



cPcyi.M=3^ j det/x(e) 
x(e) = -y j m'Ux{Q 



Cv[x(e)] 



x(e)a[x(e)] 



. (02(6, e';z) 



(02x(z) 



a[x(eO 



(A.16) 



(A. 17) 



Here we introduced the angular average (O2X of the 2D correlation (O2, associated with the 
filter Ux: 

«2x(z) = j cl0id02 Ux{Qi)Ux{Q2) (02(01, e2;z). (A.18) 

Thus, one obtains in this way the generating function cpx(j) which yields in turn the PDF 
Tx{X) from Eq.(lA31). 
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A second simple model is the "stellar model" introduced in (|290[) where we only keep the 
stellar diagrams in Eq. dA.SI) . Thus, the p— point connected correlation of the density 
field can now be written as: 

^p(xi,...,Xp) = ^ fn^2(x,,Xy). (A.19) 

The advantage of the stellar-model (|A.19I) is that it leads to very simple calculations in 
Fourier space. Indeed, Eq. (|A.19l) reads in Fourier space: 

(6(ki) . . . 5(k^)), = ^ {mfhoi^x + . . . +kp) £ Y[P{kj). (A.20) 



Using the standard exponential representation of the Dirac distribution gives (|290|) : 

as 



{XP)c=~Spf^ '^dxwP I SUx{Q)Ix{X,Qr~\ (A.21) 



where we introduced: 



(271) 

Then, using Eq. (IA.4l) we obtain: 



HX,^ = / 7^ Wx(i±1>0-) nk±,z). (A.22) 



= fdx/de t'.(e)^ <Ps (.*^;.) . (a.23) 

which directly gives cpx {y) in terms of the cumulant generating function of the 3D density 
field (p§. 



A. 3 Halo models 



An alternative to the hierarchical models presented in Section lAT2l is prov ided by the hal o 



model where the matter density field is described as a collection of halos (11761 : 12271 : 12601) . 
Then, the density correlation functions are obtained through a convolution over the halo 
density profiles. One also needs to specify the many -body correlations of the halos them- 
selves as well as their multiplicity function. Thus, the density field is written as the super- 
position of the halo profiles: 

p(x) =^m,-Mm,.(x-x/) with /dxM,„(x) = l, (A.24) 

where we introduced the normalized density profile Um{x) of halos of mass m. Then, the 
density two-point correlation function is expressed as the sum of correlations within a sin- 
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gle halo (1-halo term, both points are within the same halo) and between different halos 
(2-halo term, the two points are within two different halos): 



nil 
P 



^2(xi,X2) = J dmn{m) j dx'M„(xi - x')m,„(x2 - x') + j dmin{mi] 

dm2n(m2)^ydxiM„i(xi -Xj) Jdx2U,n2{x2-X2)^,h{A,^2^mi,m2) (A.25) 



where n{m) is the halo mass function and ^,h is the halo two-point correlation (|274l : l275h . 
In a similar fashion one can write the n-point density correlation as a sum of 1-halo up to 
/7-halo terms. In practice, it is more convenient to work in Fourier space (to simplify the 
convolution products and to take advantage of the st atistical homog eneity and isotropy of 
the system). Thus, the density power- spectrum reads (I257l : ll74l : l247|) : 



P{k)=I^{k,k) + [llik)fPt^{k) 



(A.26) 



with: 



ll{ki,...,k^)= j dmn{m)b{mf(^ u,n{h) . . .u,n{k^), (A.ll) 

where Ui„{k) is the Fourier transform of the halo density profile. In Eq. (|A.26l ) we assumed 
that the halo-halo power- spectrum can be written as Ph{k\m\,m2) = b{m\)b{m2)PL{k) 
where P^ is the linear matter power- spectrum and b{m) is the bias parameter which de- 
scribes how the halo distribution is biased with respect to the dark matter density field. 
Indeed, the 2-halo term dominates at large scales which are described by the quasi-linear 
theory (hence it is sufficient to use Pt) whereas the 1-halo term dominates a small non- 
linear scales. This also ensures that one recovers the results of standard linear theory at 
large scales. To complete the halo model one needs to specify the halo densitj^ profile, 
which is often taken from (|206|) . the halo mass fun ction , taken for instance from (i26ll:l216h. 
and the biasing of the halo distribution as from ([l99 ^. Then, the correlation functions of 
weak-lensing observables can be derived from Eqs. (|A.2l )- (IA.3l) . For instance, the 1-halo 
contribution to the angular two-point correlation of the convergence field reads (|274h : 



(K(0i)K(e2))l/,= /dx 



dm nim) 



Al^\u^{k)\'jQ{lQ), (A.28) 



where k = l/'D{x). Similar expressions give the 1-halo contribution to higher-order con- 
vergence correlation functions whereas 2-halo terms involve the linear density power- 
spectrum, 3-halo terms involve the density bispectrum and so on. This approach has been 
used to estimate the low-order corr elations of weak-lensing observables up to fourth-order 
(the kurtosis of the shear, see (l273h ) but no attempt has been made to predict the full PDF 
Tx{X) yet (however see (.27^) for a model for the tails of the matter density PDF). 
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